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In a recent letter [Phys. Rev. Lett. 108, 172003 (2012), arXiv: 1109. 2480 we have reported on 
a lattice QCD calculation of the heavy-hadron axial couplings g\, <?2, and g$. These quantities 
are low-energy constants of heavy-hadron chiral perturbation theory (HH^PT) and are related to 
the B*Bir, Ej £(, iv, and £,,* A& it couplings. In the following, we discuss important details of the 
calculation and give further results. To determine the axial couplings, we explicitly match the 
matrix elements of the axial current in QCD with the corresponding matrix elements in HrLyPT. 
We construct the ratios of correlation functions used to calculate the matrix elements in lattice 
QCD, and study the contributions from excited states. We present the complete numerical results 
and discuss the data analysis in depth. In particular, we demonstrate the convergence of SU(4\2) 
HTPyPT for the axial-current matrix elements at pion masses up to about 400 MeV and show the 
impact of the nonanalytic loop contributions. Finally, we present additional predictions for strong 
and radiative decay widths of charm and bottom baryons. 

PACS numbers: 12.38. Gc, 12.39.Fe, 12.39.Hg, 14.20. Mr 



I. INTRODUCTION 



Two important symmetries which allow many predictions in nonperturbative QCD are chiral symmetry pQ and 
Qh heavy-quark symmetry [2J. Both symmetries are a consequence of the large separation of scales in the quark masses, 
relative to the intrinsic scale of QCD. In chiral perturbation theory, an expansion is performed around the limit where 
the light-quark masses vanish, and the dynamics is determined by derivatively coupled pions, associated with the 
spontaneous breaking of chiral symmetry. In heavy-quark effective theory, an expansion is performed around the 
static limit, where a quark mass is infinitely heavy. In that limit, the heavy quark acts as a pointlike color source 
with a fixed velocity, and the spin of the heavy quark decouples. Corrections to the static limit are suppressed by 
powers of Agc^/mg, where itiq is the heavy-quark mass. Because Aq^^/to;, is particularly small, the static limit 
If-, is a good approximation to describe the dynamics of hadrons containing a bottom quark. For charmed hadrons, the 

£*-) static limit is still a reasonable first approximation. 

At the hadronic level, the combination of chiral symmetry and heavy-quark symmetry into a single effective theory 
j«~» leads to a framework known as heavy-hadron chiral perturbation theory (HH^PT), which describes the interactions 

psi of heavy-light hadrons with pions and kaons [3HH]. At leading order, the HH^PT Lagrangian contains three axial 

i—i couplings <7i , <?2 and g% . The coupling g\ determines the strength of the interaction between heavy- light mesons and 
pions, while gi and g% similarly determine the interaction of heavy-light baryons with pions. 

In the meson sector, the strong decay B* — > B n is kinematically forbidden. However, virtual pion loops contribute 
to much of the physics of B mesons, and the coupling 171 appears in calculations of important observables such as B 
meson masses, decay constants, bag parameters and form factors within chiral perturbation theory [9141 1|. Precise 
knowledge of these hadronic observables is needed for flavor physics, both within and beyond the Standard Model. In 
this context, chiral perturbation theory is needed, in particular, to describe the quark-mass dependence of lattice QCD 
results for such observables. Most current lattice QCD calculations are performed at unphysically large values of the 
up- and down quark masses to reduce the amount of computer time needed. The results from a range of quark masses 
must then be extrapolated to the physical values of the light-quark masses. Chiral perturbation theory predicts the 
functional form needed for this extrapolation, but the uncertainty in the final result is influenced significantly by the 
uncertainty in the value of the axial coupling g\ [12] . 

While most of heavy-flavor physics has traditionally focused on the B mesons, measurements at the LHC and super- 
B factories will extend the knowledge of flavor physics in the bottom-baryon sector, which provides complementary 
constraints on new physics models because of the different spin of the baryons. Therefore, lattice QCD calculations 
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of bottom-baryon observables such as form factors are needed, and as in the meson sector, chiral extrapolations of 
the data need to be performed. For baryons, the accuracy of such extrapolations can be improved dramatically if the 
values of the couplings g 2 and g 3 are known. The coupling g 3 is related to the strong decays T, y b ' — > A b ir, which are 
kinematically allowed. The widths of these decays have recently been measured at Fermilab |13) . but the experimental 
uncertainty is still large. 

The axial couplings g\, g 2l and 53 are calculable from the underlying theory, QCD. The only reliable approach 
for these nonperturbative observables is lattice QCD. While there are no previous lattice calculations of g 2 and 173, 
a number of groups have performed lattice computations of the coupling g±, both in the quenched approximation 
(i.e. neglecting the vacuum-polarization effects of light quarks) [14HT6] and with nf = 2 dynamical flavors [T7HI9] . 
In these lattice calculations, one computes matrix elements of the axial current, and relates these matrix elements 
to the coupling g\ . To fit the data and extract g\ , theoretical knowledge of the light-quark mass dependence of the 
axial-current matrix elements is required. The previous lattice calculations used fits that were linear or quadratic in 
the pion mass, or including logarithms, but the correct chiral perturbation theory predictions were not known. 

We have recently derived the expressions for the axial-current matrix elements at next-to-leading order in partially 
quenched heavy-hadron chiral perturbation theory, both for the mesons and the baryons [2U] ■ We have then performed 
the first complete lattice QCD calculation of the three axial couplings t/i, g 2 , 53, controlling all systematic uncertainties 
[2T] , In the following, we discuss important details of the analysis that were omitted for brevity in Rcf. :2l], and 
present some additional results. 

Our calculation includes nj = 2 + 1 flavors of dynamical light quarks, and makes use of data at six different values 
of the quark masses corresponding to (valence) pion masses as low as 227 MeV. Two different lattice spacings of 
a = 0.112 fm and a = 0.085 fm are used to perform a continuum extrapolation. The spatial volume is (2.7 fm) 3 , 
large enough so that finite-size effects are very small and can be removed by using finite-volume heavy-hadron chiral 
perturbation theory in the p regime. Because the axial couplings g%, g 2 , 33 are defined in the static limit, we use the 
static lattice action of Ref. [22], modified using smeared gauge links to reduce noise [23 , for the heavy quark. We 
implement the light quarks with a domain- wall action J24H26] . This is a five-dimensional formulation that realizes a 
lattice chiral symmetry for the four-dimensional theory, which becomes exact, even at finite lattice spacing, when the 
extent of the auxiliary fifth dimension is taken to infinity. As a result, the renormalization of operators is simplified 
and discretization errors are small. Our calculations make use of gauge field ensembles generated by the RBC and 
UKQCD collaborations [27] , 



This paper is organized as follows: We begin with an introduction to HH%PT in Sec. II A We derive the matrix 
elements of the axial Noether current at leading order in HHxPT using canonical quantization in Sec. |HB| The 
ratios of correlation functions used to calculate the corresponding matrix elements in lattice QCD are constructed in 



Sec. II C where we also show their spectral decomposition. The details of the lattice actions and parameters are given 



in Sec. Ill We present the numerical results for the ratios of correlation functions in Sec. IV A[ and explain our method 



for extracting the ground-state contributions to the matrix elements in Sec. |IVB| In Sec. IV C| we then describe the 
chiral fits using SU(A\2) HHxPT, including the effects of finite volume and nonzero lattice spacing. We compare our 
lattice QCD results for the axial couplings to various estimates reported in the literature in Sec. [Vj The calculations 
of heavy-baryon decay widths are presented in Sec. |VI| and we conclude in Sec. |VH| Appendix |A| contains further 
plots of numerical data, and Appendix [B] contains a comparison of our ratio method with an alternative approach 
(the summation method). 



II. AXIAL COUPLINGS FROM RATIOS OF CORRELATION FUNCTIONS 

A. Heavy-hadron chiral perturbation theory 

For hadrons containing a heavy quark (or antiquark), in the static limit tuq —¥ 00, the spin of the light degrees 
of freedom, si, becomes a conserved quantum number. The lowest-lying static-light mesons have si — 1/2, and are 
therefore grouped into pseudoscalar mesons, described by a field P l , and vector mesons, described by a field P*\ We 
work with heavy-light mesons containing a light quark of flavor i — u,d and a heavy antiquark. The vector meson field 
satisfies v^P* 1 — 0, where v is the four- velocity; this is a parameter of the effective theory, subject to the constraint 
v 2 = 1. Because of heavy-quark spin symmetry, the pseudoscalar and vector mesons are degenerate. To make the 
heavy-quark symmetry manifest, the pseudoscalar and vector meson fields can be combined into a single field H, 
which is 4 x 4-matrix- valued and given by [3] [9] 

W=[-P^+P»^} 1 —^. (1) 



This field satisfies the constraint H l -^- = H\ Next, we consider static-light baryons containing two light quarks 
of flavors i and j and a heavy quark. We include both the states with s; = and s; = 1. The states with s; = 1 
form two multiplets with J = 1/2 and J = 3/2, and are described by a Dirac spinor field B^ and a Rarita-Schwinger 
spinor field B* lJ . These s; = 1 fields are symmetric in flavor: B % i = _B Jl and B* 1 -? = -B* J \ For two quark flavors one 
has (using the notation for bottom baryons) 

B -(&^> (2) 

and similarly for B* Again, because of heavy-quark spin symmetry, the J = 1/2 and J = 3/2 baryons with si = 1 
are degenerate, and the corresponding fields can be combined into a single field [7J [55] 

s? = s# = ^(7 M + tvhsB« + b;« (3) 

satisfying U^Sff = and —^-S l J = SV. The s/ = baryons have J = 1/2 and can be described by a Dirac spinor 
field T y , which is antisymmetric in i and j and satisfies the constraint —J-T l: > = T lJ . For two quark flavors, one has 

1 ( A 6 



^VJI-A. o ')■ (4) 

In SU{2) chiral perturbation theory, the pions are described by an SU (2)-vahied field E = exp(2i$//), which 
transforms under global SU(2)l x SU(2)r transformations as 

E^LEi? f . (5) 

For the purposes of heavy-hadron chiral perturbation theory, it is convenient to also introduce the field £ = exp(z$//), 
so that E = £ 2 . The field £ transforms as 

£,^L£,U^ = U£,R\ (6) 

where the transformation matrix U(x) is a function of L, R and $(x), implicitly defined through the above equations. 
Under the vector subgroup L = R = V, the field £ transforms as £ — ► V^V"''. Therefore, the natural transformation 
laws for the heavy-hadron fields also involve the matrix U: 

s* -► u\u\s k J, 

T ij -> U\U\T kl . (7) 

The leading-order heavy-hadron chiral perturbation theory Lagrangian is then given by [5 8 



C = %-(0"Et) 5 s*< _ i tr D [B> • &H l ] - iS^v ■ $>Si j + iT ijV ■ &T* 



P 

+V2 ff3 [5^(^) 1 i T^ + Tkii^y^] + (m, terms) + (l/m Q terms), (8) 

where trrj denotes the trace in Dirac space, and the covariant derivatives are defined as 

with the vector and axial-vector fields 

fit»=± $&'£-£&'&). (10) 



The H and T fields are rescaled such that their masses do not appear in the Lagrangian. The quantity A is the 
mass difference between the S and T baryons. This mass difference does not vanish in the chiral limit nor in the 
heavy-quark limit. From experiment, one has A ra 200 MeV [131 [H]. 

Our definitions of the axial couplings gi , g 2 , and g$ in Eq. ([8]) are related to the definitions of Yan et al. [5] and 
Cho as follows [30]: 



.9i = 0i 



(Cho) 



9 



(Yan) 



„ _ „(Cho)_ 3 (Yan) 

.92 - - .92 - 2#i 



(Cho) 
.93 = .93 = 



-V3fff an) . 



(11) 



We introduced the minus sign on g 2 relative to the definition by Cho, so that in our conventions all three couplings 
are positive. 

The Lagrangian pj has the same form for both SU(2) and SU(3) chiral perturbation theory, the only difference 
being that the flavor indices run from 1 to 2 and 1 to 3, respectively. The theory can be generalized to the partially 
quenched SU(4\2) or SU (6\3) cases, where the valence and sea quarks can have different masses; for more details see 
[2"0] and the references therein. 

As can be seen by expanding the field £ = cxp(i$//) in Eq. ^ in terms of the pion field $, at lowest order the 
term with the axial coupling g\ leads to an H-H-$> vertex, the term with the coupling g 2 leads to an S-S-<& vertex, 
and finally the term with the coupling 173 leads to an S-T-<& vertex. 



B. Axial-current matrix elements in heavy-hadron chiral perturbation theory 



The simplest quantities that depend on the axial couplings are the matrix elements of the axial currents between 
single-hadron states (an alternative approach to determine the axial couplings based on static hadron-hadron potentials 
is discussed in Ref. [31 ). To extract g±, g 2 , and 53, we will calculate the matrix elements of the axial current in both 
%PT and (lattice) QCD, and match the two with each other: 



<*l 4 QCD) I^)qcd 

In QCD, the axial current is simply given by 



(X\ 4* PT > \Y) 



X PT- 



WQCD) = - < 



-fi -<* 2 "7a.759- 



(12) 



(13) 



To derive the expression for the axial current in heavy-hadron chiral perturbation theory, one can use the Noether 
procedure. For an infinitesimal local axial transformation, 



R{x) = L\x) = 1 + ia a {x)T a , 
one finds that the change in the leading-order Lagrangian is given by 

5C = (d»a a ) A°(*pt) 
with 



(14) 
(15) 



WxPT) 



*/ 2 



[r a £t^£ + r a (^£)£t _ r a m^ - r o (0 M £t) E ]« 



tvtrup^Tfjy^] 



V ^ij( T {-) k&v + Sij( T £-l k$v 

+ . gi tr D [H l {T c i + r ]1 ^w]- l g 2 e^ aX s v k y(Ti + y ] {s x y k 

+ V2g 3 [(S,)^ + r^ k + T kl {^ + r 3 {S,y k } , 
where we have introduced the quantities 

if- = \ (£V£ - £r<^) , 



TijiTf-YkT* +T ij {Tf_y h T 



■ik 



(16) 



(17) 



Equation (16) is the leading-order axial current in the chiral effective theory. In the following, we work with a 



particular flavor of the axial current, 



^-(XPT) = ^l(xPT) _ i A 2( X VT) 



(18) 



which corresponds to the QCD current d"/^ 75M. To lowest order in the pion fields (zero pion fields), the part of the 
axial current that will contribute to the matrix elements we will consider reads 

^-(xPT) = gitlD [H dl ^H u ]-ig 2 eW x v„ [(S u ) dd (S x ) du + (S v ) du (S x )™] + V2g 3 l-S^T^ - T du (S») uu } . (19) 



We will now calculate matrix elements of (191 at leading order. To this end, we canonically quantize heavy-hadron 



chiral perturbation theory. The following derivation allows us to determine the correct normalization of the matrix 
elements (canonical quantization cannot be performed in the partially quenched theory, but the normalizations can 
be inferred). We begin with the heavy mesons. Using Eq. (fTl), we find that the free part of the heavy-meson kinetic 
term in Eq. (l8|) is equal to 



C H = -iv„ tr D [H^W] 

= -iv u \-2P}d^P l + 2P* 14 8" P* 



(20) 



From this, we deduce the canonical equal-time commutation relations for the field operators (in the following, we use 
sans-serif font for operators) 



[P*(x,t),pt(x' > t)] = ^o^ 3 (x-x'), 

[P;/(x,£), PfV,*)] = -^^.«-^^)5 3 (x-x'). 

The field operators of the noninteracting theory can be written as 

d 3 fc 



1 r c\ A h 

P f (^ - — ?— f d3fc a (P)t fk)c ifc - a 

P«W~/ ( g,Ea^ ) (k, S ) £M ( S ) e --, 



where ko = v • k, and the basis polarization vectors satisfy the spin sum 

3 



E £ m( s K( s ) = ~9iu> + v H v w 



(21) 



(22) 



(23) 



Equation (21 ) is satisfied if the commutation relations for the creation and annihilation operators are 

[a* (P) (k), af )t (k')] = (2 7 r) 3 ^. ( 5 3 (k-k'), 
[a^CM), af * }t (k', s ')] - (27r) 3 ^/ ss ^ 3 (k-k'). 
We define single-particle states via 

|PXk)) = v^aP(k)|0), 
|P«(k,,s))=y2^ar )t (k, S )|0). 



(24) 



(25) 



Note that all the heavy-hadron states and operators also depend on the velocity v, which is a parameter in the 
Lagrangian. However, since v is fixed throughout this paper, we do not include a label v explicitly. 



The states ( 25 1 are normalized as 

(P l (k)|P->(k')) = 2v°{2Tr) 3 6 lJ 5 3 (k- k'), 
(P«(k,s)|P* J (k',s')) = 2 U °(2 7 r) 3 ^' ( 5 ss ^ 3 (k-k'). 



(26) 



We now calculate the matrix element (P* d \A fJ/ (x pT )|P u ). Expressing the mesonic part of the current (191 using the 
P and P* fields, we find 



^-(* PT ) = 9l tr D [H d Y To H u ] 



ffi tr D 



Pti-f-fTB (~P U K) 



1- 



*Mt nil 



= -2. 9l P7'P 



(27) 



where we have only shown the piece that contributes to the matrix element considered here. By inserting the field 
operators (22 1 into Eq. (27), we immediately obtain, at zero residual momentum, 



(P* d (0, S )|A;^ PT )(0)|P"(0))| LO = -2 5l£ ;( s ). 
Next, we consider the si = 1 baryon field S' J ' , for which the free part of the kinetic term in Eq. (tsl) is equal to 

C s = % \-iv^ + A] S'J 



(28) 



(29) 



In the following we always assume that the flavor indices on the fields S 1 ^ and Sij are ordered as i > j. We find the 
canonical anticommutation relations 



{S"«(x,t)„, S ukl p(x',t)} = - 



The field operators can be written as 



W^s^ 5 ^ 



i 



(g^-v»v v )5 3 (x-x'). 



(30) 



a/3 



s '' 0, -^7?OT/(^S a ' s » (k,s)c, " ( * K **'' 



where feo = v ■ k + 4j A, and the basis spinors satisfy the spin sum 



Y d U»{s)U v fi {s) = -{g^-v»v») 



2 



(31) 



(32) 



a/3 



Note that U^(s) is noi a Rarita-Schwinger spinor, but rather contains the degrees of freedom of both spin-1/2 and 
spin-3/2 baryons. The creation and annihilation operators satisfy the anticommutation relations 



{af s) (k )S ), aif )t (k' )S ')} = (27r) 3 <J i fc <JV. S '<5 3 (k-k'). 
Here we define single-particle states via 



|^(k, a )) = V^aif t (k,*)|0), 

which corresponds to the following normalization: 

(S ij (\c,s)\S kl (\c', s')) = v°(2Tr) 3 5 ik 6 jl 6 sa >d 3 (k-k'). 



(33) 



(34) 



(35) 



matrix element of the axial current: 

(S dd (0, S )\A^ PT H0)\S du (0, S '))\ w 



Using the expression (19) for the axial current, and the mode decomposition ([31]), we find the following leading-order 

-^g2V X e x ^U v {s)U p {s'). (36) 



Finally, we consider the si = baryon field T u . The free part of the kinetic term in Eq. &8k is equal to 






(37) 



»>j 



For the T and T fields, we also assume in the following that the flavor indices are ordered (i > j). We obtain the 
canonical anticommutation relations 



{T^(x,t) Q , T H/3 (x',i)} = J^5\8 j f 1 + r 



2v° 



k u I 



<5 3 (x-x') 



(38) 



<X0 



and the field operators can be written as 






(39) 



where k = v • k, and the basis spinors satisfy the spin sum 



Y^U a {s)U {s) = 



(40) 



a0 



The creation and annihilation operators satisfy the anticommutation relations 

{af T) (k, S ), a ^(k', S ')} = (2 7 r) 3 8 i k 8 j l S 8s ,S 3 (k-k'). 
Again, we define single-particle states via 

|T«(k, S )) = V^a|f )t (k, 5 )|0), 



(41) 



(42) 



(43) 



which are normalized as 

(T^(k,s)|T fc/ (k',s')> = v°(2Tr) 3 6 ik 6 jl 5 ss/ 5 3 (k-k'). 
Now we have all the ingredients to obtain the leading-order S-T transition matrix element of the axial current: 

(S dd (0,s)\A^ PT \o)\T du (0, s '))\ LO = -g 3 U^ S )U(s'). (44) 



To go beyond leading order, we need to replace g\, 52, 53 in Eqs. (28), (36), (44) by "effective axial couplings" (gi) s, 
(#2) off, and (53) off: 

(P* d (0,.s)|A^ PT >(0)|P"(0)) = -2 ( 5l )cff £*"(*), 
(S dd (0,s)\A^ PT H0)\S du (0,s')) = --L (fla)-ff Vxe x ^U v {a)U p {a'), 



^/2 
(S dd (0,s)\A^ PT \0)\T du (0, s ')) = -(9zUU\s)U{ S '). 



(45) 



(val) 



(val) 



(sea) 



The next-to-leading-order expressions in the partially quenched 5C/(4|2) theory (for m„ ; = rn y J al ' and rriu' 
m d a )' calculated via the perturbative expansion of the path integral, have been derived in Ref. j20j and are given 



by 



i \ ^9 X ti (vs)\ i 9\ 

(ffi)eff =9\- -pr H™>\ ') + — 



i \ 2^2 T/ ( vs ) s g 2 

(92)cB =92- -FT T \ m -n ' ) + 75 



4H(mi vs ),0)-4 4 s ^(^ vv) ,0) 
^(m( vs ),0)-4s^'(mi vv) ,0) 



+ (analytic terms), 



52 si 



/ 2 



2H(mi vs \ -A)-H(mi vv ), -A) - 2/C(m£ vs) , -A, 0) 



(analytic terms), 



(33jeff = 33 - -72" I(^Jr j ) + "rf 



H(mi vs \ -A)--H(ml vv \ -A) 



+\H{m^\ A) + 3H(mi vs ), A) - JC{m^\ A, 0) 



33 g\ 



U{m%«\ A) - T^m^, A) +H(m£ s \ 0) - ^^(m^, 0) 



+ (analytic terms). 



(46) 

Here, mi denotes the mass of a pion consisting of a valence and a sea quark, mi denotes the mass of a pion 
consisting of two valence quarks, and 5y S ~ [mi vv ^] 2 — [mi ] 2 . The functions X, H, H^ and /C, which arise from 
the chiral loops and include the leading effects of the finite spatial volume, are defined in Ref. [2D]- At the order 
considered here, the analytic terms in Eq. (461 are linear functions of [mi vv) ] 2 and [mi ] 2 . 



C. Axial-current matrix elements in lattice QCD 



For the lattice QCD calculation, we construct interpolating fields for the heavy hadrons in terms of the quark fields 
as follows: 

P l =Qaa (75) a/9 £,3, 
P ? = Qaa (Tittup <7a/3> 



V" e abc (C 7m )^ 7 q a/3 q bl Q ca 
T % J = e abc (C7 5 )^ 7 q l a p q J bj Q ca 



(47) 



Here, a, 6, c are color indices, and a, /3, 7 are spinor indices. The light-quark field of flavor i is denoted by q l (we will 
also use the notation u = q u , d — q d ), and the heavy-quark (antiquark) field is denoted by Q (Q). The tilde on the 
light-quark fields indicates that these are smeared over multiple spatial lattice sites, in order to improve the overlap 
of the interpolating fields with the corresponding ground-state hadrons and reduce excited-state contamination. We 
use gauge-invariant Gaussian smearing obtained by 



2 \ ™s 

— A< 2 ) 



(48) 



where A' 2 ' is a three-dimensional gauge-covariant lattice Laplacian which includes gauge links, a is the smearing 
width and ns is the number of smearing iterations. 

The heavy quark Q is defined in the static limit, and we set v = 0. The static heavy-quark field Q satisfies 



1 + 7 C 



Q = Q. 



(49) 



Note that the interpolating field S$ a couples to both the J = 1/2 and the J = 3/2 baryons with si = 1, with exactly 
the same relative amplitudes as |3|. 

We use the local four-dimensional lattice axial current, where the quark and antiquark fields are evaluated at the 
same lattice point. This current requires a finite renormalization factor Za to match the continuum QCD current: 



^a(LQCD) = Za g ioQ l(r°)< .( 7M75 ) a/3 g^. 



(50) 



The value of Za depends on the lattice action used, as well as the lattice spacing. We will come back to this in 



Sec. Ill As before, we will work with a specific flavor of the axial current: 

^-(LQCD) = Z A d aa (^ l5 ) a/3 U a0 . 



(51) 



In the following, we will omit the superscript "— (LQCD)" in the axial current (51 1. Next, we define the overlap 
factors, which describe the overlap of the operators corresponding to the interpolating fields (47) with the relevant 



ground-state hadrons in QCD. Here, we use the same notation and normalization of states as in Sec. II B 

(0|P*(0)|P*(fc))=Z pi , 

(0\P; i (0)\P* i (k,s))=Z P , i e„(s), 

(0\S%(0)\S ij Oc,a))=Z 3 i i U lia {a) ) 

(0|ri(0)|^(k, S ))=^ T «W a ( S ). (52) 

We stress that these states are now meant to be hadron states in (lattice) QCD, rather than in the chiral effective 
theory. 



We calculate Euclidean two-point functions of the interpolating fields (47), as well as Euclidean three-point functions 



with an insertion of the axial current (51 1 . These Euclidean correlation functions are obtained from the lattice path 



integral, which is performed numerically using importance sampling. In the following, we assume that the Wick 
rotation t — > —it has been performed, so that t denotes the Euclidean time. 

We calculate the following three-point functions, where (...) denotes the path-integral over the gauge- and fermion 
fields (for the domain-wall action used in this work, there is also an additional path integral over Pauli-Villars fields 



c\p* d a Pir^ o = EE < p * dM ( x -*) a v &J) ^(o) >, 

x x' 

c[s m a s du rj(t, o-EE< s^(x,t) ^(x',o s p duP (o) ), 

x x' 

C[S dd A T du ]^{t, t')=J2J2< ^"(x,i) ^(x',0 T dufj (0) ), 

x x' 

C[T d ^ S dd ]%(t, t')-EE( T * l ( x ^)^ t ( x '^')^ /3 (0)). 



(53) 



x x' 



In addition, we calculate the two-point functions 

ClP u Pt](t)=J2(P u (x,t)Pt(0)), 

X 

c[p* d P^T u {t) = E ( P* df2 (*, t) p* "t(o) ), 

X 

C[S dd s dd C (t) = J2 ( s?"(x,t) s dd0 (o) >, 

X 

c\s d - s du Q(t) = e ( sf(x, t) r duP (o) ), 

X 

c[T du T du u(t) = J2 ( r Q du (x,i) T duP (0) ). 



(54) 



As an example, we show how the two-point function C[S •3 dd ] l ^ l/ *(t) is constructed in terms of quark propagators. 
The lattice calculation is performed in the isospin limit of equal up and down quark masses. Inserting the definitions 
of the interpolating fields, and performing the Grassmann integrals over the fermion fields explicitly, we have 



-~a , -~S 



C[S dd S M ]%(t) = e abc (Cy% e fgh {Cr) P „ E ( #(*) %(*) <£(*) d p (0) d„(0) Q p (0) 

X ^ 

= e abc (Cy% e fgh (Ur) pa J2 ( <^>, °) <% )w (*, 0) Gf Q)a ,(x, 0) 

X 

- %)^(^0) <^ 7 >,0) G$ )ap (x,0)) D 
= 2 e abc (CY% e fgh (CTU E (&&sM °) &%)„(*> 0) G$ )afi (x, 0) 



(55) 
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Here, Gig) denotes a smeared-source smeared-sink light-quark propagator and G(q) a heavy-quark propagator. The 
notation ( ... )u indicates the path integral over the gauge fields U only. In the last step in Eq. (55) we have used 



the symmetry of (Cjfi) and the antisymmetry of e a bc to equate the two Wick contractions. Note that the two-point 
function C[S du Sdu]^g(t) contains only one Wick contraction, because the two light quarks have different flavors. 
Therefore, in the isospin limit, 



For the static lattice action we are using, the heavy-quark propagator is equal to 



G{Q) a p{x,Q) = <5x, 



1 + 7° 



a0 



t^(x,i - a) Ufat - 2a) ... t/ o f (x,0), 



(56) 



(57) 



where Uq is a spatially smeared gauge link in the temporal direction (more details will be given in Sec. III). Because the 
static heavy-quark propagator (57) contains the Kronecker delta <5 x ,o> the sums over x in all the correlation functions, 
Eqs. (53) and (54), are eliminated. To calculate the three-point functions (|53|) in terms of quark propagators, we use 



pairs of light-quark propagators with smeared sources at (0, 0) and (0, t). This means that new inversions are needed 
for each value of t. 

By using the Hamiltonian and the momentum operator to shift the left interpolating operator from x to and the 
axial current from x' to 0, and inserting complete sets of states, we can show that 



1 3 

C[P* d APlr(t, t') = -Z P , d Z P ^ e*(«) (P* d (0, S )|A"(0)|P"(0)> 



e -E p , d (t-t') e -Eput' 



s=l 
6 6 



C[S dd A S du f a 7(t, t') = Z SM Z* d „ J2 E U ^ s ) (S d \0,sW(0)\S du (0 lS i)} U P p {s') e 



E s dd(t—t ) „~E sdu t 



s=l s' = l 
6 2 



C[S dd A TtoCpfa t') = Z sdd Z* Tdu J2J2 U ^ (S dd (0,sW(0)\T du (0,s / ))U^s / )e- E s^ t ~^ e' E ^ 1 ' + 

s=l s' = l 
2 6 

c[T du A i s dd fj;(t, t') = Z T ^Z* sdd ^2^2u a (s) (T du (0,s)\A^(0)\S dd (0, S ')) U" (s') e - £ W*-0 e -E sdd t' + 



s=l s' = l 



(58) 



and 



C[P™Pt){t) = -\Z P u?e 



2 „-E Pu t 
3 



e -E p , d t 



C[P* d P*Jr(t) = l -\Z P , d \^e^s)e*%s) 

s=l 

C[S dd S dd ]^(t) = J2 U ^ TjV ^\ Z S' 

s=l 

C[S du S du ]^(t) =J2 U ^ T7 ^ l^-l 2 e~ Esd 

s=l 

2 

C[T du T du } afs (t) = Y^U^U^s) \Z TdM 



2 e -E sdd t + ___ ^ 



2 p — E Tdu t 



(59) 



where the ellipsis indicates contributions from excited states, whose contribution relative to the ground-state contri- 
bution shown here vanishes exponentially for t — >• oo, if — > oo, i — t' — > oo (here and in the following we assume an 
infinite temporal extent of the lattice). We will return to the discussion of excited states at the end of this section. 



Using the relations ( 23 ) , ( 32 ) , and ( 40 1 to perform the spin sums (recall that here we have v = 0) and Eq. ( 45 1 to 
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express the axial-current matrix elements in terms of the effective axial couplings {gi) c s, (g%) e &, (<?3) e ff, we obtain 
C[P* d APlf v {t, t') = ~ ( ffl ) off («V - cf u ) Z P , d Z* pu e - B p*4t-*') e - Eput ' + ... , 



C[S dd A S du C/(t, = -^ ( 52 ) cff «a 



z \yuvp 



1 + 7° 



Z S ddZ% iu e' B s dd( - t - t "> e - E s^ t ' + 



C[S M A T du ]^(t, t') = -(<7 3 )eff (vV - <T) 

C[T d Mt S rfrf Qi, t') = -( 53 ) cff (vV - fT) 



a/3 

1+7 01 



2 
1+7 01 



a/3 



a£ 



Z gddZrpdu 6 S 



E s dd{t—t') p — E T d u t 



e~"Td«" + .. 



Zj'duZ g dd e T ' 



E T d u {t—t ) „ — E s ddt 



e -^sdd>. _l_ __ 



(60) 



and 



C [ P *d P*t]^(£) = -(vV - 5 ^) |Z P , d | 2 e-Ep*^ + 



'1+7° 



[1 + 7°] 


2 

"1 + 7°" 


2 



\Z S dd\ 2 e -^ ddt + ... , 



a/3 



|Z<j,^ e"-^"* + 



C[T " Tci u ] a( g(i) 



a/3 



(61) 



a/3 



In the following, we remove the trivial spin-structure [ "y ° 1 _, which comes purely from the heavy-quark propagator 



(571, from all baryon correlation functions. 



Because the lattice calculation is performed in the isospin limit (and in the static limit for the heavy quark), we 
have the relations 

Epu = Epd* , 

Zpu = Z pd* , 

Egdd = Egdu , 

Z Qdd = v2 Zgdu 



(62) 



[the factor of -\/2 in the last line comes from Eq. (56)]. As a consequence of the equality of energies, the i'-dependence 
of the ground-state contribution in the three-point functions C[P* d AP^"(t : t') and C[S dd A S du ]^ p (t, t') cancels 
completely. For these three-point functions, the ^-dependence as well as the Z factors in the ground-state contribution 
can be canceled by forming the ratios 



Ri(t, t') 



\Y,l=iC[P* d APlY»{t,t') 



C[P u Pl]{t) 



(63) 



(.9 



lJcff 



and 



o (+ ^ o b K,u, P =i ^ P c ^ ddA s du r p (t, f) 

H 2 (t, t ) = I 



lEUc[s^s dd ]^(t) 



(64) 



(ff2)off 



where, as before, the ellipsis indicates contributions from excited states that vanish exponentially when all Euclidean 
time separations are taken to infinity. To extract (g3) c ff, we use the double ratio 



R 3 (t, f) = 



1 v^3 



I EjUi C[S dd A T du ]w(t, V) | Y.%,1 C[T<*Mt S M ]w(t, f) 



\ 



(53) 



1 v^3 



E^=i C\S dd S dd ]^(t) [C[T d « T du ](t)] 



3 ^m 



(65) 



cir 
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The numerical results for (63), (64), and (65) and the subsequent analysis will be described in Sec. IV In the following, 



we discuss the contributions from excited states to the ratios. Again, we assume an infinite temporal extent of the 
lattice; with a finite temporal extent T this means that the following discussion is only valid for source-sink separations 
t that are smaller than T/2 by a sufficient distance (which is the case in our numerical calculations). We begin with 
R\, and define 



(0|P«(0)|P«) = Z F , n , 

K d (0)\P: d (e)) = Z 



P,n Efii 



CCX(0)|^(e)H-2A| 



pp*) 



(66) 



where \P%) denotes the n-th excited state with a nonzero overlap (0|P u (0)|P^), and similarly for \P* d (e)). Because 
of heavy-quark symmetry and isospin symmetry, all energies and Z-factors in the P u sector are equal to those in the 



P* d sector, and A {P S ] = A*^"'. Note that 



4?^ = (SiW 



The complete spectral decomposition of R\ reads 



Y^oo y^°° 7 7* /t( PP *) -E P „i „-E P m (t-t') 

D f . ,,s l^ n =\l^ m =\Zp.;nZ Pm A nm e p -" e F -^ > 

l( ' } ~ Y-°° \7 I 2 -E P t ' 

Showing only the contributions from the ground states and first excited states, we find that 

■M 2 ( A (P p ') _ a {pp ' ] \ p -'p* + z P . lZp2 ( P p.) Spt , z Pa z* P1 .*(pp*) s P (t-t') 

, A | 1^22 -^11 I c ' \Z p ,i\ 2 12 e ' \Z P1 \ 2 ^12 e 



(67) 



(68) 



Ri(t, t') = A[ pp 



-Zp,2 

Zf 



\ Zp, 2 \ 2 -8 P t 

I Zp ; i I 



(69) 

with the energy gap Sp = Epp — Ep s \, For a given value of t, the smallest contamination from excited states is 
obtained at the midpoint t' = t/2. Evaluating (69) at t' = t/2, we get 



R 1 {t,t/2)=A[ PP ' ) 



\Z Pj2l \1 , A (PP*) _ APP*)^ -S P t 
\Z P A ^22 -^11 



V 



23? 



Zp.iZ p2 .(PP*) 

\z P A 2 12 



e -iSpt 



1 



*P,2 Vjpj 



Z p ,: 



where 3? denotes the real part. By using the Taylor expansion 1/(1 , 7 

2 

^22 



g-op 

Z P , 2 I 2 -<5 P i 



e 



) = 1- 



Zp t2 \ 2 „-8 P t 



Ri(t, t/2) 



.(pp*) 



ni 



Z 



P. 2 



z 



p,l 



(4 PP * 



A 



(PP*)\„-8 P t 



11 



)e" 



23? 



■Zp,l^p,2 .(pp*) 
I^P.ll 2 ^ 



e -J«P* 



(70) 

we obtain 

(71) 



where we have omitted terms that decay like e 2°p t or faster, and are therefore exponentially suppressed relative 
to the terms shown in (71) at large t. The result for R2(t, t/2) has the same form (with suitable definitions of the 



overlap factors and matrix elements): 



R 2 (t, t/2) = A[i s) + 



Zs.2 

Zs.i 



(A 



(SS) 
22 



A (SS) )e - Ss t 



23? 



z s,iZ* S2 
llZs^F 



a(ss) 

A \2 



■iSst, 



(72) 



with A\ x ' — (<?2)cff and 8s = Eg 2 — Eg.i- For the double ratio Rs(t, t/2) we obtain, after Taylor-expanding the 
square root and omitting terms that decay faster than e~ Sst or e~ Srrt , 



Sa(M /2) = A<f»-' 4 ' f>1 ^ 1 ^' ^T^-*' 



+3? 

+3? 



2|Z s ,i| 2 21^12 


AST) Zs,2 

A 2\ ry 

Z S\1_ 


erhsst + % 


AST) Z T ,2 
A \2 7 

Z/ T,l_ 


e^ s 


AST) Z T,2< 

22 ry 


7* 

°S,2 


e -h(Ss + 


8 T )t + 





(73) 



with A 



(ST) 



(33)cff and 63 = Eg 2 — Eg t \, St = £p,2 — E-ri- The ^-dependent terms in the right-hand-sides of Eqs. 



(71), (72), and (73) are the leading excited-state contributions to the extraction of (<7i)eff from the ratio method. 
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III. LATTICE ACTIONS AND PARAMETERS 



Our calculations are based on gauge field ensembles generated by the RBC/UKQCD collaboration. These ensem- 
bles include 2 + 1 dynamical light-quark flavors, implemented with a domain- wall action [24Tf26] , The gluons are 
implemented with the Iwasaki action [321 I33j , which is known to reduce the residual chiral symmetry breaking of the 
domain- wall action [53] . 

The details of the actions used in generating the ensembles can be found in Ref. [35] . Our analysis includes the 
ensembles of size 24 3 x 64 and 32 3 x 64, which are described in Ref. [27]. These ensembles have lattice spacings of 
approximately 0.112 and 0.085 fm, respectively, so that the spatial volume is about (2.7 fm) 3 in both cases. 

At the coarse lattice spacing, we used only the ensemble with arny^T — 0.005, which is the lightest available 
mass. At the fine lattice spacing, we used the ensembles with the lightest two available values of the sea quark mass, 



(sea) 



0.004 and am^ s J ) 



0.006. The values for the residual quark mass, which is the additive quark-mass 



am 

renormalization coming from the residual chiral symmetry breaking at finite extent of the fifth dimension L s , are 
approximately am ros = 0.0032 at the coarse lattice spacing and am Tes = 0.00067 at the fine lattice spacing. The 
sea-strange-quark masses are about 10% above the physical value [27], and we assign a 1.5% systematic uncertainty 
to our final results for the axial couplings to account for this, based on the size of the effect on similar axial-current 
observables as determined using mass reweighting in Ref. [27j . 

We calculated light-quark propagators using exactly the same domain-wall action that was used by the RBC/UKQCD 
collaboration for the sea quarks, with the same domain- wall height of aM§ = 1.8 and extent of the fifth dimension 
L s = 16. We used propagator sources smeared according to Eq. (48), with a = 4.35 and rig = 30. As summarized 
in Table I we calculated propagators for valence quark masses omj^ both equal to and lighter than the sea quark 

masses. The data with m^ & d < m^ are referred to as "partially quenched" . Also shown in Table I and plotted 
in Fig. [I] are the corresponding masses of pions composed of the three different possible combinations or valence and 
sea quarks. The lightest valence- valence pion mass is 227(3) MeV, at the fine lattice spacing. 



i 3 xr 



(scaj 



(sea) 



(val) 



a (fm) 



mi ss) (MeV) 



,W 



(MeV) 



,(w) 



(MeV) 



24 3 x 64 


0.04 


0.005 


0.001 


0.1119(17) 


336(5) 


24 3 x 64 


0.04 


0.005 


0.002 


0.1119(17) 


336(5) 


24 3 x 64 


0.04 


0.005 


0.005 


0.1119(17) 


336(5) 


32 3 x 64 


0.03 


0.004 


0.002 


0.0849(12) 


295(4) 


32 3 x 64 


0.03 


0.004 


0.004 


0.0849(12) 


295(4) 


32 3 x 64 


0.03 


0.006 


0.006 


0.0848(17) 


352(7) 



294(5) 
304(5) 
336(5) 

263(4) 
295(4) 
352(7) 



245(4) 
270(4) 
336(5) 

227(3) 
295(4) 
352(7) 



TABLE I. Parameters of the gauge field ensembles and quark propagators. The lattice spacing values are from Ref. |36| . 



> 

CD 

a 



0.15 



0.10 



S 0.05 



0.112 fm 
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[GeV 2 ] 



0.15 



,(w)i2 



FIG. 1. Fhe values of [raj- ] and [m4 vs ] used in our calculation. The dashed line indicates the unquenched case rrvk 



(vv) 



,(vs) 



For the heavy quark, we use a static action of the form given by Eichten and Hill [22 , which corresponds to heavy- 
quark propagators of the form (57). For the temporal gauge links in this action (or, equivalently, the propagators), 
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we use HYP smeared links [37] with smearing parameters ot\ — «2 = «3 = 0.75. This leads to an exponential 
improvement in the signal-to- noise ratio |23j . The smearing procedure can be iterated tjhyp times, leading to a 
broader smearing and further improvement of the signal-to-noise ratio. We generated data for «hyp = 1, 2, 3, 5, 10. 
While all of these actions have the same continuum limit, the dependence of the results on the lattice spacing is 
expected to be different for different values of hhyp- One may naively expect large discretization effects for large 
values of tihyp > which correspond to a large spatial extent of the heavy-quark-gluon interaction vertex. We will discuss 
this in detail when giving the numerical results in Sec. |IV| Our final axial coupling results only make use of data 
generated with Uhyp = 1, 2,3. 



As mentioned in Sec. II C to perform the contractions for the three-point functions (53), we required pairs of light- 
quark propagators with sources located at the same spatial point and separated by tja steps in the time direction. 
The numbers of measurements (propagator pairs) for each value of t/a are given in Table III] At the coarse lattice 
spacing, our data come from typically 120 statistically independent gauge field configurations; at the fine lattice 
spacing we used about 240 statistically independent gauge configurations of the arajj c ^ = 0.004 ensemble and 180 for 

the amf^JJ" = 0.006 ensemble. In most cases, we have more measurements than configurations, because we generated 
propagators from multiple source points on the lattice. In those cases, we have averaged over the source locations 
before the further analysis to remove possible autocorrelations. 



L 3 xT 


(val) 
am u,d 


t/a 






iVmeas (apprOX.) 


24 3 x 64 


0.001 


10 






550 


24 3 x 64 


0.001 


9,8, 


7, 


6 


240 


24 3 x 64 


0.001 


5 






460 


24 3 x 64 


0.001 


4 






120 


24 3 x 64 


0.002 


10 






880 


24 3 x 64 


0.002 


9,8, 


7. 


6,4 


240 


24 3 x 64 


0.002 


5 






480 


24 3 x 64 


0.005 


10 






960 


24 3 x 64 


0.005 


9,8, 


7. 


6,4 


240 


24 3 x 64 


0.005 


5 






480 


32 3 x 64 


0.002 


12 






1200 


32 3 x 64 


0.002 


9, 6 






480 


32 3 x 64 


0.004 


12 






1200 


32 3 x 64 


0.004 


9, 6 






480 


32 3 x 64 


0.006 


13 






700 



TABLE II. Number of propagator pairs used for the three-point functions for various values of the source-sink separation t/a. 



Within each of the three gauge field ensembles that we used, the data from different source-sink separations, 
different valence quark masses, and different values of tihyp & r e correlated with each other. In our analysis, we 
properly took into account these correlations using the statistical bootstrap procedure. The initial averaging over 
source locations mentioned above was also required to reduce all data from the same ensemble to matching ordered 
sets of measurements, as necessary to calculate the covariance matrices. It turned out that the correlations between 
the am d — 0.001 and nm^' = 0.002/0.005 data at the coarse lattice spacing were very weak even though the 
data came from the same ensemble of gauge field configurations. The reason was that all source locations used for 
the am^j = 0.001 propagators were distinct from those used at amj^' = 0.002/0.005. In contrast, the data from 

am ud ~ 0.002 and am'j = 0.005 came from almost identical source locations, resulting in very strong correlations 
(these correlations were advantageous in constraining the quark-mass dependence in our chiral fits). Similarly, at the 
fine lattice spacing, the data from am J = 0.002 and arri^d ~ 0-004 came from identical source locations, leading 
to strong correlations. 

For the axial-current renormalization parameter, we use the values obtained nonperturbatively by the RBC/UKQCD 



collaboration, which are 
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Z, 



0.7019(26) for a = 0.112 fm, 
0.7396(17) for a = 0.085 fm. 



(74) 



IV. DATA ANALYSIS 



Ratios 



Examples of numerical results for the ratios (63), (64), and (65) arc shown in Fig.[2](for a = 0.112 fm) and Fig. [3] 
(for a = 0.085 fm). These ratios were calculated using statistical bootstrap to take into account the strong correlations 
between the three-point and two-point functions in numerator and denominator. Because of these correlations, the 
statistical uncertainties in the ratios are found to be smaller than those in the three-point functions themselves. To 
maximize correlations, it is essential to use the two-point functions from the source locations matching those of the 
three-point functions. The figures show results at two different values of the source-sink separation t in each case, for 



am 



(val) 



0.002 and tt-hyp = 3. 
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FIG. 2. Ratios Ri, R%, R3 as a function of the current insertion time slice t', for the source-sink separations t/a ■ 
t/a = 10 (right). The data shown are for a = 0.112 fm, tihyp = 3, and awi„™ = 0.002. 



5 (left) and 



We observe clear plateaus of Ri(t, t') as a function of t' , whose extent appears to be slightly larger for the simple 
ratios R\ and R2 compared to the double ratio R3. The plateaus in if indicate that in these regions contributions from 
off-diagonal matrix elements of the axial current between a ground-state hadron and an excited state are negligible in 
comparison to the statistical uncertainties, because such a transition matrix element would introduce ^'-dependence 
[see Eqs. (68) and (69)]. For R\, the flatness with respect to if has previously been observed and discussed in 

Refs. [mnsr 

We averaged Ri(t, t') over a symmetric range of if values around tjl in a region where there was no statistically 
significant if dependence (requiring that the x 2 /d.o.f . of correlated constant fits be of order 1). These regions and the 
extracted values, which we denote as Ri(t), are indicated in Figs. [2] and [3] for representative data sets. The averaging 
in the plateau region is essentially equivalent to using 



Ri(t, t/2). 



(75) 



Indeed, because of the strong corrrelation between neighboring if points, we found that the averaging over if in the 
plateau region (which is again performed using bootstrap) gives almost the same result and uncertainty as Ri(t, i/2). 
An alternative method for defining Ri(t) is discussed in Appendix |B| 

Summary plots of all the extracted values for i?i(i), i?2(i), and R3[t) at the different quark masses, lattice spacings, 
and numbers of heavy-quark smearing iterations iijjyp are given in the appendix in Fig. |18| Note that the numbers 
of measurements vary (see Tabic Ujh. The statistical uncertainties are found to grow quickly when t is increased or 
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FIG. 3. Ratios R\, R2, R3 as a function of the current insertion time slice if , for the source-sink separations t/a — 9 (left) and 
t/a = 12 (right). The data shown here are for a — 0.085 fm and ?ihyp = 3, and (urn^d — 0.002. 



,(val) • 



"u,d 



is decreased, as expected 



Furthermore, the statistical uncertainties are reduced with every iteration of 



HYP smearing in the static heavy-quark action, which is also expected [23]. While the results for (gi) e & — lim^oo itj(i) 
from all fixed values of iihyp will become equal when the lattice spacing is taken to zero, at nonzero lattice spacing 
different values of tihyp will have different discretization errors. Indeed, R\ and R2 are seen to increase with riHYP- 
Remarkably, the results for R3 are almost independent of 71hyp within the statistical uncertainties even at nonzero 
lattice spacing. See also Fig. [6]m Sec. IVB for the nnYP-dependence of (gi) c e = lim^oo Ri(t). 



B. Extrapolation to infinite source-sink separation 

The effective axial couplings (<?i) e ff at given values of the lattice spacing, the quark masses, and iihyPi are defined 
as the infinite-time limits of Ri (t) : 



(ft) 



cir 



lim R t (t). 



(76) 



The i-dependence of Ri(t) is caused by excited states, and at sufficiently large t, the contributions from the first 
relevant excitation dominates. The expected functional form of i?i(i), R2(t), and Ra(t) in this regime is shown in 
Eqs. (71 1, (72 1, and (73), respectively. As can be seen in these equations, the "diagonal" contributions proportional 



to the matrix elements An and A22 decay exponentially like 



(77) 



where S is the energy gap to the first excited state that has a nonzero overlap with the corresponding interpolating 
field as defined in Eq. (47); 5 = 5 P for i? 1; S = 63 for R 2 , and S = mm(5s, S T ) for R 3 . Additionally, there are 



"off-diagonal" contributions proportional to the matrix elements A12 and A21, which decay only like 



(78) 



However, as discussed in Sec. |IV A and Ref. [18 , these off-diagonal contributions appear to be small, because the 
numerical results for R(t, t') show plateaus as a function of if . 

We performed correlated fits to the lattice data for Ri(t, a, m, «hyp), which depends on the source-sink separation 
t, the lattice spacing a, the quark masses (here denoted generically as m), and rtHYP, using the following functions: 



i?i(i, a, m, n HY p) - (ffi)rf(o, m, n HY p) - A^a, to, n HY p) e -*t w,!W »' 1 
R 2 (t, a, to, n HY p) = (.9 2 )cff (a, m, n HY p) - A 2 {a, to, n HY p) e - s ^ a ' m > n ^ *, 
R 3 (t, a, to, n HY p) = (fl3)rfE(o, to, n HY p) - A 3 (a, to, n HY p) e ~ 5 ^ m > n ^ *. 



(79) 
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Here Si is the energy gap to the dominant excitation in Ri. Because Si is not constrained, it does not matter whether 
or not we include a factor of 1/2 in the exponent. Since the energy gap Si is positive by definition, we choose to 
parametrize it as 



a 5i(a,m, riHYp) 



_ Ji(a,m,n HY p) 

f i 



(80) 



using the logarithm Zj(a,77i,riHYp) as the fit parameter. 

Because the statistical uncertainties in Ri grow exponentially as t is increased, we were only able to perform the 
lattice QCD calculations in the range t < 1.1 fm. As can be seen in Table fill at the coarse lattice spacing (a = 0.112 
fm) we have data for t/a = 4,5,6,7,8,9, 10. We found that the functions (79) described the data from the coarse 
lattice spacing well for all these values of t/a. The smallest statistical uncertainties are obtained for am™ d = 0.005 
and riHYP = 10, and therefore we first performed unconstrained fits to this data set, obtaining the following fit results 



(li converted to Si — 

(3l)cff 
(32)off 
(53)eff 



,h 



/a): 

0.5264 ±0.0090, A x 

1.037 ±0.033, A 2 

0.827 ±0.032, A. 



0.53 ±0.39, 
0.73 ±0.35, 
0.98 ±0.25, 



Si 

S2 

s 3 



(1.08 ±0.38) GeV, X^/d.o.i. 
(0.75 ±0.29) GeV, vVd.o.f. 
(0.66 ±0.16) GeV, Y 2 /d.o.f. 



1.07, 
0.31, 
0.41. 



(81) 



The corresponding fits and the data points are shown in the left-hand side of Fig. [4J Notice that the fit to R\ gives an 
energy gap consistent with recent lattice QCD results for the 2S radial excitation energy of about 0.75 GeV [39|l40], 
This is expected, because our interpolating fields (47) for the heavy-light mesons have S- wave-type smearing and 



therefore do not couple to the lower-lying IP state. The fit result (and the flatness of the plateau as a function of t') 
indicates that the data for R\ are dominated by 

2 



i?i(i)«(ffi)rfE + 



J P(2S) 



J P(1S) 



(PP*) 

with a negligible off-diagonal matrix element A^ s 1S ~ 



( A (PP') 
\^-2S,2S 



,(PP' 



A 1S,1S ) e 



-Set 



(82) 
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FIG. 4. Fits of the ^-dependence of Ri(t), R-2(t), Rs(t), for tihyp = 10, using the functions defined in Eq. (791. The left-hand 



side shows unconstrained fits of the a — 0.112 fm, am 



(val) 



0.005 data; the right-hand side shows fits of the a = 0.085 fm, 



fval) 



0.004 data, where the energy gap was constrained using information from (81 1 



At the fine lattice spacing, we only have data for t/a = 6,9,12 (for am^ d = 0.002, 0.004) or t/a — 13 (for 
rf Ktd = 0-006). Because the energy gap in physical units is not expected to have an a-dependence that is larger 



than the statistical uncertainties in (81), it is possible to use the fit results for U from the coarse lattice spacing to 



constrain the parameters li at the fine lattice spacing (at similar pion masses). As a first step, we performed fits to the 
data from the fine lattice spacing with amj^ = 0.004 and «hyp = 10, where only the parameters li were constrained 
using an augmented \ 2 with a Gaussian prior for /,, 



{li - h? 



(83) 
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Here, U and ui i are the central values and uncertainties of the (scaled) energy gap parameters from the fit to the coarse 

0.004 



am u,d = 0.005, tt-hyp 
and riHYP = 10 then gave 



(val) 



10 data, Eq. (81). The fits to the data from the fine lattice spacing with am u d 



(ffi)«ff = 0.510 ±0.011, Ai 
(<? 2 )eff = 0.998 ±0.041, A 2 

(g 3 )cff = 0.805 ± 0.034, A 3 



2.7 ±3.6, *i 

1.8 ±1.7, S 2 
2.0 ±1.0, S 3 



(1.08 ±0.38) GeV, 
(0.75 ±0.29) GeV, 
(0.66 ±0.16) GeV. 



(84) 



At the fine lattice spacing, it was necessary to remove the data points with the shortest separation t/a — 6 to obtain 
acceptable single-exponential fits. Therefore, the resulting gap parameter matches exactly the prior, and ^ /d.o.f. 
is undefined. The central values of the overlap parameters Ai in (84 1 are larger than in (81), indicating a stronger 



overlap of the interpolating fields with excited states at the fine lattice spacing. Different overlap factors were expected 
here, because the smearing width of the light-quark fields in physical units was different (we used the same smearing 
width in lattice units for both lattice spacings). The fit curves are shown in the right-hand side of Fig. El 

We then performed new fits to the data for all values of the quark masses and njjYP- For those fits, the parameters 
(5i)cff (a, to, riHYp) were left unconstrained, but Gaussian priors were used for both li and A4, with central values and 
widths as taken from the initial fits (81) for the coarse lattice spacing and (84) for the fine lattice spacing. Examples 
of these fits are shown in Fig. [5J The only assumption made by using the priors is that the dependence of li and Ai on 
"■hyp and on the quark masses is smaller than the width of the priors as determined by the statistical uncertainties 
in (81 ) and (|84|. Given that these widths were 25% or larger in all cases, this appears to be a reasonable assumption. 



As a test, we also performed unconstrained fits where pos sible , wh ich gave consistent results but were less stable. 
The results for (<7i) e ff(o, »to, ^hyp) are gi ven i n Tables III and IV and plotted in Fig. pi The central values and 
uncertainties shown there are bootstrap averages and 68% widths, respectively, from repeated correlated fits of the 
t-dependence for a bootstrap ensemble of data. 
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FIG. 5. Fits of the t-dependence of Ri(t), R,2(t), R-i(t), for jihyp = 3, using the functions defined in Eq. 

0.002 data, and the right-hand side shows the a = 0.085 fm, arn^ 



(val) 
u,d 



1. p}. 

(vaTj— _ 



side shows the a = 0.112 fm, am 

both cases, the overlap parameters Ai and energy gaps 6% were constrained using information from (81 1 and (841 



The left-hand 
0.002 data. In 



To estimate the systematic uncertainties caused by higher excited states, we calculated the shifts in (<7i) e ff for a 
representative data set at the coarse lattice spacing (omj,™' = 0.002, ?ihyp = 3) when removing one or two data 
points with the smallest t/a (= 4, 5) from the fits, or adding a second exponential to the fit function, 



Ri(t) = (j9i)«B - A-i e 



-Si t 



B, A; 



-(*i+*f ) )t 



(85) 



(2) 

Because the available data were not sufficient to determine the new parameters B t and S^ , we used Gaussian priors 
to constrain these parameters to physically reasonable values. The parameters Bi describe the amplitudes of the 
second-excited-state contribution, relative to the first-excited-state contribution, and we set Bi = 0, e rg. = 2. For the 



9(2 



energy gaps a8\ = e * we used priors with central values equal to 2/3 times the fit results for <5.; in (81 1, and widths 



of 100%. The fitted parameters (si)cff an d the corresponding shifts 5(gi) e g for the different cases are"shown in Table 
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(sea) 
am u,d 


(val) 
am u,d 


«HYP 


(Sl)cff 


(ff2)eff 


(ff3)eff 


0.005 


0.001 


1 


0.463(28) 


1.14(16) 




0.005 


0.001 


2 


0.473(20) 


1.094(92] 




0.005 


0.001 


3 


0.479(18) 


1.077(74; 


0.843(68) 


0.005 


0.001 


5 


0.488(15) 


1.063(62; 


0.822(50) 


0.005 


0.001 


10 


0.514(15) 


1.075(54; 


0.828(40) 


0.005 


0.002 


1 


0.499(17) 


0.984(42; 


0.815(51) 


0.005 


0.002 


2 


0.496(13) 


0.996(35; 


0.816(41) 


0.005 


0.002 


3 


0.499(11) 


0.993(29; 


0.810(37) 


0.005 


0.002 


5 


0.5059(96) 


1.001(27; 


0.814(35) 


0.005 


0.002 


10 


0.5230(89) 


1.039(35; 


0.828(36) 


0.005 


0.005 


1 


0.496(13) 


0.986(36; 


0.831(43) 


0.005 


0.005 


2 


0.4950(94) 


0.986(25; 


0.820(34) 


0.005 


0.005 


3 


0.4987(80) 


0.990(23; 


0.812(30) 


0.005 


0.005 


5 


0.5080(74) 


1.006(2i; 


0.814(28) 


0.005 


0.005 


10 


0.5270(71) 


1.039(25; 


0.828(27) 



TABLE III. Effective axial couplings (gi) c a at a — 0.112 fm, obtained by extrapolating Ri(t) to t — Vt 

we do not have results for jihyp = 1,2, because the statistical fluctuations were too large to calculate the square root of the 



oo. At am^' = 0.001 



double ratio, Eq. ( 65 I 



(sea) 
l u,d 



(val) 
b u,d 



"HYP (ffl)eff (ff2)eff 



(.93 



0.004 


0.002 


1 


0.496(52; 


0.95(11) 


0.78(14) 


0.004 


0.002 


2 


0.507(3i; 


0.920(75) 


0.788(85) 


0.004 


0.002 


3 


0.505(24; 


0.940(66) 


0.778(64) 


0.004 


0.002 


5 


0.501(17; 


0.946(59) 


0.762(49) 


0.004 


0.002 


10 


0.505(13; 


0.962(50) 


0.772(41) 


0.004 


0.004 


1 


0.488(38; 


0.939(82) 


0.799(85) 


0.004 


0.004 


2 


0.498(23; 


0.948(65) 


0.795(62) 


0.004 


0.004 


3 


0.502(18; 


0.982(66) 


0.803(49) 


0.004 


0.004 


5 


0.503(14; 


0.995(53) 


0.799(39) 


0.004 


0.004 


10 


o.5ii(io; 


1.001(41) 


0.807(33) 


0.006 


0.006 


1 


0.412(49; 


0.86(16) 


0.79(12) 


0.006 


0.006 


2 


0.452(33; 


0.905(95) 


0.807(73) 


0.006 


0.006 


3 


0.465(29; 


0.925(81) 


0.797(64) 


0.006 


0.006 


5 


0.481(26; 


0.974(72) 


0.805(60) 


0.006 


0.006 


10 


0.508(23; 


1.030(68) 


0.824(56) 



TABLE IV. Effective axial couplings (gi) c a at a = 0.085 fm, obtained by extrapolating Ri(t) to t — oo. 



\V\ Since the shifts 5(gi) e s themselves have statistical uncertainties crS(gi) e g, we choose to quote the maximum value 
of 



V[ S (g*)cs} 2 + W5(9i)cs] 



(86) 



from the three different methods (removing t/a = 4, removing t/a = 4, 5, or adding a second exponential) as our 
estimate of the systematic uncertainties in (ffi)eff caused by higher excited states. The final estimates are 1.7%, 2.8%, 
and 4.9% for g\, <?2, and 173, respectively. 

Because the shifts 5(gi) e g in Table M are consistent with zero in most cases, an alternative way of estimating the 
systematics is to consider only the increase (calculated with quadrature) in the uncertainties of the fitted parameters 
when the fits are modified by removing data points or including higher-order terms. We will use that method for the 
chiral fits at the end of Sec. IV C see Eq. ([96]) . The two different methods (86) and (96 1 for calculating the size of 
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FIG. 6. Plot of the results for 



G fr from Tables 



III 



and 



IV 



For each combination of a and am J*. , results from up to five 



different values of «hyp are shown (from left to right: tihyp = 1, 2, 3, 5, 10; points offset horizontally for legibility) 



the systematic uncertainties give consistent values. 



Fit 


(ffl)cff 


<%l)eff 


(fl^cff 


<%2)cff 


(93) off 


S(gs) e g 


Original 


0.499(11) 





0.993(29) 





0.810(36) 





t/a — 4 removed 


0.496(13) 


0.0030(76) 


0.975(35) 


0.016(19) 


0.783(43) 


0.026(15) 


t/a — 4, 5 removed 


0.494(12) 


0.0041(76) 


0.984(41) 


0.009(26) 


0.807(54) 


0.003(30) 


Second exponential added 


0.498(11) 


0.0009(77) 


0.988(30) 


0.005(21) 


0.796(40) 


0.014(36) 



TABLE V. Fits used to estimate systematic uncertainties from higher excited states. Data for a = 0.112 fm, am u 



(val) 



0.002, 



hhyp = 3. Shown are the fit results for (<?i) e ff, as well as the differences 5(gi)eB to the original fit result, calculated using 
bootstrap. 
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C. Extraction of the axial couplings g±, 32, and gz using HHyPT fits of the data 



In the previous section, we obtained results for the effective axial couplings {gi) e s{a,m, 71hyp) at two different 
lattice spacings a, multiple values for the quark masses ara^™ and amj^ , and multiple values for the heavy-quark 
gauge- link smearing parameter njjYP (corresponding to multiple heavy-quark lattice discretizations) . All data are for 
a finite spatial volume of about (2.7 fm) 3 . In the following, we discuss how we extracted the axial couplings g\, g 2 , $3, 
which are the parameters of the continuum heavy- hadron chiral perturbation theory Lagrangian ([8]), from the data 
for (g i ) cS (a,m,n iiY p). 

To fit the quark-mass and volume dependence of (gi) c s we use the next-to-leading-order predictions from HHyPT 
[20] . which were already shown in Eq. (46). Here, we extend these formulas to include the leading effects of the 
nonzero lattice spacing a. These leading effects are quadratic in a, with coefficients di inilYP that depend on tihyp- 
We do not expect O(a) errors because of the chiral symmetry of the domain- wall action used for the light quarks 
(neglecting the small effects caused by the residual chiral symmetry breaking at finite L s [57]) and the automatic 
0(a) improvement of the static heavy-quark action. Higher-order effects in the a- and m^-dependence of (gi) e s are 
discussed at the end of this section. 

For (gi)eff, the fit function is given by 



(ffi)«ff(0) m, n H yp) = 9i 



1 



^I(mH) + ^4H( m ^,0) 



. c (w) [m (vv) ]2 + c (v S ) [m (v S ) ] 2 + di] 



■4<$ s ft, 



'{™ ( r\ o)} 



2 
"HYP a 



(87) 



Similarly, for (g 2 ) e s and (g 3 ) e s, we use 



(32)eff(a, m, TJ-HYP) = .92 



1 - 



l(m^) + |f {- H(m£ s \ 0) - 5 2 vs H vl {m^\ 0)} 



+ ^{2-H(m£ B \ -A)-H(m<T\ -A)-2/C(m 



+ C 2 VV) [m^? + 4 VS) [m { : s) ] 2 + d 2 , „ HYP a 2 



(vs) 



0)} 



(53)eff (a, m, n H Yp) = .93 



Arn^) + ^{H(m^\ -A) - \%{mt\ -A) 
+ \u{m^\ A) + ZU{m^\ A) - IC(m^\ A, 0)} 



P 
3, 



+ 



<J 2 



{-H(m™, A)-H(m£ v \ A) + n(m£ s \ 0) - 6^ s H v >(m^, 0)} 



(w) [ ra H,2 , >) r TO (vs)i2 



+ 4 VV) K 



^3, riHYP ^ 



(89) 



The functions I, H, T-L,y and K. are the chiral loop corrections [20 . They include finite- volume effects and therefore 
they also depend on the lattice size. Furthermore, these functions depend on the renormalization scale fx, but this 



„(w) 



scale-dependence is absorbed by the fit parameters c\ and c\ , as we checked explicitly by varying /j, in the fits. 



(vs) 



We set the pion decay constant to / = 132 MeV and the S - T mass splitting in Eqs. ( 88 ) and ( 89 1 to A = 200 MeV 



This value of A is consistent with experiments [T51 12"5] and with our lattice data (we also checked that varying A 
within a few percent does not significantly affect the results for the axial couplings). We calculated the covariances of 
all correlated data points in Tables [Hi] and [TV] using bootstrap, and performed fully correlated fits using the inverse 
of the covariance matrix in the definition of x z . This method propagates the uncertainties and correlations of {gi) c s , 
as obtained from the fits to the ratios Ri, into the extracted parameters gi of the HHyPT Lagrangian. 

Results from fits of the (gi) fi data using the function (87) are given in Table VI and Fig. [71 The fit parameters are 



9i, c 



(vv) (vs) 



and {d 



1, "HYP 



} (the latter for all values of 71hyp that were included in the fit). We performed fits that 
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included data with multiple values of "hyp, as well as individual fits including only data with one value of "hyp- The 
fits that included multiple values of «hyp enforced a common continuum limit of the data with different "hyp , but 
with separate a 2 -coefficients d\, „ HYP for each ii H yp- While we know that the actual continuum limit for all values of 
"hyp has to be the same (if we took a to zero in the numerical calculations) , we only have data for two different values 
of a, and one may question whether the approach of the continuum limit is described by a simple a 2 -dependence as 
assumed in Eq. (87 1. In particular, one may be worried that large values of "hyp, which correspond to more spatially 
extended heavy-quark actions, could lead to non-negligible contributions from higher powers of a [23] . To investigate 
this, we started from a fit that included all values of "hyp (1, 2, 3, 5, 10), and then successively removed the data 



with the largest values of «hyp ■ As can be seen in Table VI and Fig. [7J the fit including the data from all values of 
"hyp had a poor quality, Q = 0.17, and gave a somewhat low value for g\. After excluding "hyp = 10 and njjYP = 5, 
the fits had a good quality and the results for gi were stable under further exclusions of the largest 7iHYP _vam es. The 
fit including tt-hyp = 1, 2, 3, which has Q = 0.70, gave the result 



.gi =0.449 ±0.047, 



(90) 



Estimates of the systematic uncertainties in ( 90 1 will be given at the end of this section 
including only one value of "hyp were all consistent with (90), even for whyp = 10 
of a 2 are actually negligible for the values of the lattice spacings considered here 



The results from the fits 
This suggests that higher powers 
a = 0.085 fm and a = 0.112 fin). 



The deviating value of g\ as well as the poor Q for the fit that included all values of "hyp simultaneously are likely 
caused by technical issues with the covariance-fitting of highly correlated data, associated with small eigenvalues of 
the data correlation matrix [5T] . We will return to the discussion of higher-order discretization effects at the end of 
this section. 



"HYP 


ffi 


d.o.f. 


xVd-o.f. 


Q 


1, 2, 3, 5, 10 


0.371(28) 


30-8 


1.3 


0.17 


1, 2, 3, 5 


0.401(39) 


24-7 


1.2 


0.29 


1, 2, 3 


0.449(47) 


18-6 


0.75 


0.70 


1,2 


0.440(60) 


12-5 


0.85 


0.54 


10 


0.450(38) 


6-4 


0.09 


0.91 


5 


0.468(47) 


6-4 


0.61 


0.55 


3 


0.482(55) 


6-4 


0.73 


0.49 


2 


0.465(66) 


6-4 


1.0 


0.36 


1 


0.49(10) 


6-4 


0.72 


0.49 



TABLE VI. Results for the mesonic axial coupling gi, obtained by fitting the data for (<?i) ff using the function (871. The first 



four rows show the results from fits which include data with multiple values of the heavy-quark smearing parameter tihyp- 
The remaining rows show the results from fits with only one value of hhyp- The number of degrees of freedom (d.o.f.) is given 
in the form (number of data points) — (number of fit parameters). The last column of the table gives the quality of the fit 
Q = T (d.o.f./2, x 2 /2). 



For the baryonic axia l coup ling s, we performed simultaneous, fully correlated fits to the data for {g2)eS and (<?3) ff 

using the functions (88) and (89), with the fit parameters c/2, <?3, <4 , C3 , c 2 , cij and {d2,n H YP> <fe, «hypj (the 
latter for all values of "hyp that were included in the fit). As already discussed in the fits for (<7i) e ff, we performed 
fits that included data with multiple values of "hyp , as well as individual fits including only data with one value of 
"hyp- The results are shown in Table VII and Fig. [8] Again, we select the fit that includes «hyp = 1,2,3, which 
gives 



52 = 0.84 ± 0.20 s tat, 
.93 = 0.71 ± 0.12 stat . 



(91) 



This fit also had the highest value of the quality of fit, Q = 0.92. Estimates of the systematic uncertainties in (91) 



will be given at the end of this section. As can be seen in Fig.[8l the results (91 ) are in fact consistent with the results 
from all other fits within the statistical uncertainties, demonstrating that heavy-quark discretization errors are under 
good control. The covariance matrix for gi and 53 is 



Cov = 



0.040 0.011 
0.011 0.014 



(92) 
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FIG. 7. Graphical representation of the fit results for gi from Table |VI) The horizontal axis corresponds to the different fits, 
ordered (from left to right) in the same way as the rows in the table (from top to bottom). The line and shaded region in the 
upper plot indicate the selected result and its uncertainty, which is taken from the third fit (the fit that includes data with 
IlHYP = 1,2,3). 



"HYP 



<)> 



()■>, 



d.o.f. 



xVd-o.f. 



1, 2, 3, 5, 10 


0.72(12) 


0.635(90) 


58- 


- 16 


0.94 


0.57 


1, 2, 3, 5 


0.73(13) 


0.61(11) 


46- 


- 14 


1.1 


0.31 


1,2,3 


0.84(20) 


0.71(12) 


34- 


- 12 


0.61 


0.92 


1,2 


0.81(22) 


0.57(17) 


22- 


- 10 


0.50 


0.91 


10 


0.90(15) 


0.75(11) 


12- 


-8 


0.64 


0.64 


5 


0.98(19) 


0.76(13) 


12- 


-8 


0.74 


0.57 


3 


0.98(23) 


0.74(15) 


12- 


-8 


0.54 


0.71 


2 


0.91(23) 


0.66(18) 


12- 


-8 


0.51 


0.67 


1 


0.79(29) 


0.61(27) 


12- 


-8 


0.42 


0.74 



TABLE VII. Results for t he b aryo nic axial couplings <?2 and 173, obtained by simultaneously fitting the data for (<72)eff and 
(j?3)cff using the functions (88 1 and (891. The first four rows show the results from fits which include data with multiple values 
of the heavy-quark smearing parameter tihyp. The remaining rows show the results from fits with only one value of riHYP- 
The number of degrees of freedom (d.o.f.) is given in the form (number of data points) — (number of fit parameters). The last 
column of the table gives the quality of the fit Q = V (d.o.f ./2, x 2 /2) ■ 



The corresponding likelihood function is plotted in Fig. |9] 

As another check, we performed fits of (gi) c g where we excluded all the partially quenched data (i.e keeping only 
the unitary data with m™ = m4 ). In that cas e, onl y one ana lytic counterterm is needed for each coupling, and 



we removed the terms 



-O) r™( vs )i2 



from Eqs. (87 1, (88), and (89). These fits, again using tihyp = 1,2,3, then gave 



g x = 0.467 ± 0.056, g 2 = 0.92 ± 0.22, and g 3 = 0.72 ± 0.14, inftill agreement with Q and J91J) and with slightly 
larger uncertainties. 

Plots of the functions (gi) c g, (32)09, and (<73) G ff, with the parameters from the fits including the complete data 
the fit that gives the results (90) and (91 1] are shown in Figs 



with riHYP = 1,2,3 [i.e. 



10 



11 



and 12 



For the 
0.085 fm, and 



figures, the functions were evaluated in infinite volume, for the lattice spacings a — 0.112 fm, a 
a — 0. The right-hand sides of the figures show the values and uncertainties of the fitted functions for the unitary case 
mi = mJr , while the left-hand sides show the dependence on both m^ and rrvk ■ At the two nonzero values 
of a, the functions were evaluated for jihyp = 3 and the corresponding data points are also shown (in the continuum 
limit, the functions for njjYP = 1, 2, 3 are all equal). To allow the inclusion in these plots, the data points were shifted 
to infinite volume using 

(9i)eff,data("i, £ = oo) = (fl>i) e ff,data(w, L = 2.7 fm) + [(si) e ff,fit (m, L = 00) - (ft) e ff,fit(w, L = 2.7 fm)] , (93) 

where we use the notation m = (m,x , m; ) . The numerical values of the volume shifts are given in Table VIII The 
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FIG. 8. Graphical representation of the fit results for c/2 and g% from Table |VII| The horizontal axis corresponds to the different 
fits, ordered (from left to right) in the same way as the rows in the table (from top to bottom). The lines and shaded regions in 
the upper two plots indicate the selected results and their uncertainties, which are taken from the third fit (the fit that includes 
data with jihyp = 1, 2, 3). 
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FIG. 9. Likelihood function for gi and (73, equal to L( 



.,«» 



2 ,g 3 ) = (27T)" 1 det(Cov)- 1 / 2 expj-i( Si - ^[Cov" 1 ]^ - gf)} 



where g\ are the central values of our fit results (91 1 and Cov is the covariance matrix (921. The dashed curve indicates the 
standard error ellipse. 



largest volume shift (2.8%) occurred for (<72)eff at m^ = 227 MeV 



rn ( J s) (MeV) 


mi vv) (MeV) 


(9l)iff J -(w)^ J 


(92)^' 


(s 3 )ia , -(ff 3 )i^ 


294 


245 


0.0057 


0.015 


0.0074 


304 


270 


0.0040 


0.0070 


0.0027 


336 


336 


0.0016 


0.00037 


-0.00079 


263 


227 


0.0072 


0.028 


0.013 


295 


295 


0.0031 


0.00027 


-0.0012 


352 


352 


0.0013 


0.00033 


-0.00071 



TABLE VIII. Size of the finite-volume corrections for the pion masses where we have data. 
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11 



and 



The functions (<72)eff and (<73) e ff develop small imaginary parts for pion masses below the S —> Ttt threshold at 
TOjr = A [20] (the lattice data are all above the threshold). The extracted parameters 31,2,3 are real. Figures 
12 show the real parts of (<?2)eff and (<73) G ff only, which have kinks at the thresholds. 

The fit results for the parameters c\ vv , c\ , which describe the analytic contributions, were natural-sized, i.e. of 



order 1/Ar. with A x f=a inf^, for the renormalization scale /i = 4jvf v . The fit results for the parameters di 



which 



describe the lattices-spacing dependence, were also of natural size and consistent with zero within the statistical 
uncertainties. The absence of significant a-dependence can also be seen in Figs.[lOj[ll] and[L2| 

The individual contributions from different classes of Feynman diagrams in HH^PT (20) to t he fitted functions 

■' ys ') are shown in Figs. 



13 



14 



Note that 



(gi)eff, 3t[(<72)cff], and ^t[{gj,) c a] (evaluated for a = 0, L = 00, and mX = m* 
while the sum of all contributions (including the analytic terms) is independent of the renormalization scale /z, the 
individual contributions are not, and the figures are based on the natural scale \x = 47r/ ff . For the range of pion masses 
considered here, the NLO contributions are significantly smaller than the LO contribution (which is equal to gi). This, 

and the natural size of the fitted coefficients c\ vv , c^ , indicates that the chiral expansion of the axial-current matrix 
elements is well-behaved here. 

To estimate the size of systematic uncertainties caused by the missing NNLO terms in the fits to the quark-mass 
and lattice-spacing dependence, we performed fits to the data using modified functions (<?i) cff that include 

higher-order analytic terms: 



/ n(NLO+HO)/ v / s(NLO), n 



+gt 



v) [m^] 4 + Cl (vs ' vs) [m^)] 4 + Cl (vv ' vs) [m[ vv )] 2 [m( vs )] 



,(vv) 



,(w)l2 , j(vs) 2 

i ir J T "i, tihyp 



i { : s) ? + ^ 



(94) 



Here, the functions (ffj) jg are as defined in Eqs. (|87|), (|88|), and (|89|). Because we do not have enough data to fit all 



the parameters in Eq. (94), we constrained the parameters corresponding to the higher-order terms using Gaussian 



priors centered around zero and with widths equal to some dimensionlcss factor w times the relevant natural scales: 

(vv,vv) n , / . 4 

c\ ' — ± w/A 



(vs,vs) 
C i 

(vv.vs) 

j(w) 



,(vs) 



± w/A 4 , 
± w/A 4 , 
± w A 2 QCD /A 



d^ HYP =0±^A 2 QCD /A^ 

hi, mhyp — ± w Aq CD . 



(95) 
Here we used A x = Airf^ with /„. = 132 MeV, and Aqcd = 300 MeV. The fit results for the axial couplings gi as a 



function of the width factor w are given in Table |IX| While the case w = corresponds to the original NLO fits, in 
the limit w — > 00 the new parameters would become unconstrained (because we have insufficient data, we are unable 
to perform fits in this limit). 



fli 



5a(gi 



',12 



Sa(g 2 



!)-', 










0.449(47) 





0.84(20) 





0.71(12) 





1 


0.449(47) 


0.0020 


0.84(20) 


0.0023 


0.71(12) 


0.0045 


5 


0.452(48) 


0.0089 


0.84(20) 


0.014 


0.70(12) 


0.017 


10 


0.455(50) 


0.016 


0.84(20) 


0.024 


0.70(12) 


0.026 


50 


0.464(72) 


0.054 


0.82(22) 


0.099 


0.68(15) 


0.094 


100 


0.452(94) 


0.082 


0.78(26) 


0.17 


0.63(21) 


0.17 



TABLE IX. Results of higher-order fits using Eq. ( 94 1 as a function of the width factor w defined in Eq. ( 95 1 



As can be seen in Table |LX| the shifts in the central values of the axial couplings are smaller than the statistical 
errors up to the very large width w = 100. This is a consequence of the smallness of the quantities m 4 /(47r/ 7r ) 4 , 



2G 

a 2 AQ CD m^/(47r/ 7r ) 2 , and a 4 Aq CD for the pion masses and lattice spacings where we have data. The shifts in the 
central values fluctuate statistically and can be close to zero even for large w (at least for gi). However, including the 
higher-order terms leads to a systematic increase in the uncertainties of the fit parameters gi (as calculated from the 
Hessian of x 2 )> as expected. Also shown in the table is the quantity 



8a( 9i ) = ^(^(NLO+HO)-^.) 



(NLO) 



(96) 



where cr(gi)( NL °^ is the original uncertainty oi_Qi from the NLO fit, and cr(f/i)^ NLO+HO - ) is the new uncertainty of gi 
from the hig her-order fit (95). To calculate (96l we used more digits for cr(<?i) ( - NLO " ) and cr(<7i)^ NL0+H0 ^ than shown 
in Table IX Equation (961 gives the additional uncertainty in gi, calculated using quadrature, that results from the 
higher-order terms. This additional uncertainty 5a(gi) scales roughly linearly with the width parameter w. For a 
reasonable choice of w, the quantity 5a{g{) can be considered to be the systematic uncertainty in gi from the NLO 
fit due to the missing NNLO terms. Here we choose the conservative value of w — 10 for this purpose. The resulting 
estimates of relative systematic uncertainties can be found in Table [X] There, we also show the estimates of the other 
relevant sources of uncertainties: effects of higher excited states in the fits to Ri(t) as discussed in Sec. IV B and the 
effects the sea-strange-quark mass being about 10% above the physical value, as discussed in Sec. |IH| Including the 
estimates of the total systematic uncertainties, our final results for the axial couplings, based on (90 1 and (91), are 
then 



9l = 0.449 ±0.047 s tat ± 0.019 syst 
52 = 0.84 ±0.20 stat ±0.04 syst 
.93 = 0.71 ±0.12 stat ±0.04 syst 



= 0.449 ±0.051, 
= 0.84 ±0.20, 
= 0.71 ±0.13. 



(97) 



Source 



9i 



y-2 



Unphysical value of m 
Total 



(sea) 



<)', 



NNLO terms in fits of rn-n- and a-dependence 
Higher excited states in fits to Ri(t) 



3.6% 


2.8% 


3.7% 


1.7% 


2.8% 


4.9% 


1.5% 


1.5% 


1.5% 



4.2% 4.3% 6.3% 



TABLE X. Estimates of systematic uncertainties in the axial couplings 
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FIG. 10. The fitted function (gi) c s, evaluated in infinite volume, for whyp = 3, at different lattice spacings (from top to 
bottom: a — 0.112 fm, a — 0.085 fm, a — 0), along with the data points (shifted to infinite volume). The left-hand side shows 
the dependence on both mi and m; T . The right-hand side shows the function (and its statistical uncertainty) evaluated at 



,(w) 



,(vs) 



(vv) 



,(vb) 



In the plots on the right-hand side, the partially quenched data points, which have m\ < mi- , are indicated 
with open symbols. They are shown at m n = mi- , even though the fitted function (gi) B B actually has slightly different values 
for these points. 
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FIG. 11. Like Fig. 10 but for the real part of (<?2)eff- 
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FIG. 12. Like Fig. 10 but for the real part of (<?3) e ff. 
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FIG. 13. Contributions from individual classes of Feynman diagrams in HHxPT (see Ref. 20 ) to the fitted function (gi) B s, 
evaluated in the infinite volume, continuum limit, for mj v — rrVh - The renornialization scale is /i = 4irf. 
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FIG. 14. Contributions from individual classes of Feynman diagrams in HHxPT (see Ref. 20 ) to the real part of the fitted 

functions (<?2)eff (left panel) and (33)08 (right panel), evaluated in the infinite volume, continuum limit, for mi = mi . The 
renornialization scale is fj, — 4vr/. 
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COMPARISON WITH PREVIOUS RESULTS FOR THE AXIAL COUPLINGS 



We begin this section by discussing previous lattice calculations of the heavy-meson axial-current matrix elements 
and the corresponding extractions of 171 . A summary of results is shown in Table |XI| All of the past works used an 
order-a improved Wilson action [32 for the light quarks, and variants of the Eichten-Hill action [22, 23] for the static 
heavy quark. The first lattice estimate for g\ was obtained in the pioneering work of Ref. [13] , using a 12 3 x 24 lattice 
and quenched gauge fields, where the fermion determinants in the path-integral weight are set to 1, which means 
that the vacuum-polarization effects of the light quarks are neglected. In Ref. [14], the average of (gi) e s from two 
different valence pion masses (760 MeV, 900 MeV) was taken as the result for g 1 . Quenched calculations of g\ were 
also reported in Refs. [TS] and [IB]. The results for g\ in these works were obtained by extrapolating data for (gi) c s, 
at pion masses in the range of about 550 to 850 MeV, linearly in [mi ] 2 to m4 = 0. 



Reference 



n/, action 



,Hl2 



(GeV 2 



9i 



De Divitiis et al, 1998 Q3] 
Abada et al, 2004 [15] 
Negishi et al, 2007 [16] 
Ohki et al, 2008 [17] 
Becirevic et al, 2009 [18] 
Bulava et al, 2010 [19] 

This work 



0, clover 
0, clover 
0, clover 
2, clover 
2, clover 
2, clover 



0.58- 
0.30- 
0.43- 
0.24 
0.16 
0.063 



0.81 
0.71 
0.72 
1.2 
1.2 
0.49 



0.42 ±0.04 ±0.08 

0.48 ±0.03 ±0.11 

0.517 ±0.016 

0.516 ± 0.005 ± 0.033 ± 0.028 ± 0.028 

0.44 ± 0.03t°:oo 

0.51 ±0.02 



2 ±1, domain wall 0.052-0.12 0.449 ± 0.047 stat ± 0.019 syst 



TABLE XL Comparison of lattice QCD results for the mesonic axial coupling g\. Also shown are the numbers of dynamical 
light-quark flavors n/, the fermion lattice action, and the range of valence pion masses used in the calculation. 

Since calculations without sea quarks have uncontrolled systematic errors, more recent lattice calculations of g\ 
have been performed with dynamical flavors, albeit only for n/ = 2. The first of these was done in Ref. [17j . using 
two different lattices of sizes 12 3 x 24 and 16 3 x 32, and pion masses in the range 490 - 1100 MeV. Stochastic all-to-all 
propagators were used to reduce the statistical uncertainties. In Ref. |17| . the data for the axial-current matrix elements 
(ffi)cff w as fitted using different approaches: linear in m 2 , linear±quadratic in m 2 , or linear±quadratic±logarithmic 
in m 2 , using the average of g\ from the linear and the linear±quadratic±logarithmic fits as the final result. A second 
unqucnchcd calculation was published in Ref. [TB] , using three different lattice spacings and pion masses in the range 
from 400 to 1100 MeV. In Ref. [18], the coupling g\ was obtained from a linear±logarithmic fit of (<?i) c ff. Recently, 
the axial couplings of orbitally excited heavy-light mesons were also included [33] . Another n f = 2 calculation of g\ 
was reported in Ref. [ID], with three different lattice spacings and pion masses down to 250 MeV. In Ref. [TDJ, the 
result of an extrapolation of (gi) c ff linear in to 2 was given as the value of g\. 

The coefficient of the chiral logarithm used in the fits of the axial-current matrix elements (<7i) G ff in Refs. [17] 
and [IS] was set equal to that of the strong decay P* — » P it in SU(2) HHxPT [33] , because the corresponding 
loop contributions for the axial-current matrix elements were not known at that time. The NLO expression for the 
strong-decay amplitude in SU(2) HHxPT is proportional to 



M(P* ^Ptt)cx 5 i 



1- 



n (^uy 



log- 

[i 1 



(98) 



We have recently derived the NLO expressions for the axial-current matrix elements in SU(2), SU(3), SU (4\2), and 
SU(6\3) HLLyRT [2D]. As discussed in Ref. [3D], the chiral expansion of the axial-current matrix elements contains an 
additional tadpole loop contribution, which modifies the coefficient of the logarithm. In the SU(2) case, one has 



(5i) 



cir 



,9i 



1 - (2 + 4. 9 2 ) 



(^uy 



log- 



It' 



(99) 



Because g\ 

of the logarithm in Eq. 

(see Fig 



0.5, the coefficient of the logarithm in Eq. (99) is numerically about 3 times larger than the coefficient 

(98). The logarithm makes (gi) c fi as a function of ra 2 curve downward when m 2 is decreased 

10|). The results for g\ from the previous unquenched lattice calculations, which incorrectly used Eq. (98) 



or did not include any logarithm in the fits, would be significantly lower if the correct HHxPT formula (99 1 had been 
used instead. We have attempted fits to the data of [TTHTD] using Eq. (99), obtaining values of g\ that are about 10 
to 20% lower than what is published in these works. Note, however, that HHxPT is not expected to converge in the 



32 



upper range of the pion masses in [TTJ [IS] • For the data used in the present work, incorrect fits linear in m% or using 
Eq. (98) give values for g\ that are higher than the correct result, Eq. (90), by 12% and 8%, respectively. 



Next, we move to the discussion of various theoretical estimates of the axial couplings gi, g 2 , 53 based on approxi- 
mations, models, and experimental data. A comparison of these estimates to our QCD results is shown in Table |XIl) 
The nonrelativistic quark model (NRQM) predicts g\ — g^ d , <?2 = 2g^ d and <?3 = \f i 2,g v ^' [5J, where g^ = 1 is the axial 
coupling of the single-quark transition u — > d. Interestingly, if g^ 1 is set to 0.75, the value needed to reproduce the 
experimental value of the nucleon axial charge, one obtains g\ = 0.75, g 2 = 1-5, 53 = 1.06, still significantly larger than 
our QCD results. The predicted ratios of the axial couplings in the NRQM are, however, consistent with our lattice 
determination. The relativistic quark models of Refs. [35] and [IS] give g\ = 1/3 and g\ = 0.6 ± 0.1, respectively. 



Reference 




Method 


3i 
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93 


Yan ei al., 1992 [5] 




Nonrelativistic quark model 


1 


2 




V2 


Colangelo ei al, 1994 [431 




Relativistic quark model 


1/3 








Becirevic, 1999 gB] 




Quark model with Dirac eq. 


0.6 ±0.1 










Guralnik ei al, 1992 [47] 




Skyrme model 




1 


6 




1.3 


Colangelo ei al, 1994 [48] 




Sum rules 


0.15- 0.55 










Belyaev ei al, 1994 g5] 




Sum rules 


0.32 ±0.02 










Dosch and Narison, 1995 [50] 


Sum rules 


0.15 ±0.03 










Colangelo and Fazio, 1997 


m 


Sum rules 
Sum rules 


0.09 - 0.44 










Pirjol and Yan, 1997 El 


<^ 




9l 


< V2 


Zhu and Dai, 1998 [53] 




Sum rules 




1.56 ±0.30 ±0.30 


0.94 ±0.06 ±0.20 


Cho and Georgi, 1992 [54] 




B[D* ->Dir],B[D* -4D7] 


0.34 ±0.48 








Arnesen ei al, 2005 El] 


B[D\ a 


^fl (s) i],B[i)^D( s) 7],f[D*] 


0.51 








Li ei al, 2010 [56] 




dr[B -> ■ntv] 


< 0.87 








Cheng, 1997 [30] 




T[E* ->A c 7r], NRQM 


0.70 ±0.12 


1.40 ±0.24 


0.99 ±0.17 


This work 




Lattice QCD 


0.449 ±0.051 


0.84 ±0.20 


0.71 ±0.13 



TABLE XII. Comparison of our lattice QCD results for the axial couplings g\, gi, and 33 with other determinations as reported 
in the literature [all results are shown in our normalization, see Eq. ( |11[ )]. Here, NRQM stands for nonrelativistic quark model. 
Where decay widths or branching fractions are listed under "Method" , these are experimental inputs. As discussed in the 
main text, the axial couplings extracted from experimental data are defined away from the static limit in some cases. When 
a reference contained multiple results for the same coupling and did not specify which one is the most reliable, we quote here 
the range from the lowest result minus its uncertainty up to the highest result plus its uncertainty. 



Another theoretical approach for estimating the axial couplings is the large- N c limit of QCD, where N c is the 
number of colors. In the limit N c — > 00, one finds that the baryonic couplings satisfy the relation [4"Tll5"?] 



53 



Af c =oo 



1.22. 



(100) 



For comparison, our lattice QCD result for this ratio is 



(101) 



with results as shown in Table XII For the 



— = 1.19(26), 

53 

and the nonrelativistic quark model predicts 52/53 = v2 « 1.41. 

The axial couplings have also been estimated using sum rules 
heavy-meson coupling g±, most sum rule determinations are smaller than our lattice QCD result, and much smaller 
than the NRQM value. In contrast, the values of gi and 33 obtained using sum rules in Ref. [53j are larger than our 
lattice results. 

Experimental data for various heavy-hadron decay processes has also been used to determine the axial couplings. 
In Ref. [54] . electromagnetic interactions were included in HHxPT, and the coupling g\ was extracted from the 
measured branching fractions B[D* — > D ir] and B[D* — > D 7] at tree level, finding g\ — 0.43 ± 0.61 for m c = 1.5 
GeV and g\ = 0.34 ± 0.48 for m c = 1.7 GeV. Note that these values for g\ are not defined in the static limit; they are 
effective values corresponding to the D*Dir coupling. A similar calculation, which additionally included the leading 
nonanalytic effects in the radiative decays, is reported in Ref. |58j . The complete I/ttiq and loop corrections in both 



:« 



the strong and radiative decays were included in the analysis of Ref. [35] . There, the fit to experimental data for the 
branching fractions B[D1 s -, — > D^ tt] and B[D1-, — > D^ 7] gave two possible solutions for g±. The fit of Ref. [59] was 
updated later by including experimental results for T[D*] [60], leading to g\ ~ 0.51 [55], where (unlike in Ref. [59] ) 
gi is defined in the static limit. 

Recently, g\ was also extracted from data for the B — > ixlv form factors, giving results for g\ in the range from 
0.02 ± 0.32 up to 0.73lg '^2 depending on the parametrizations of the form factor shape [SS]. The measured widths of 
the baryonic decays E* — > A c 7r were used in Refs. [52] and [30] to estimate 33, with the result 0.99 ±0.17. The NR QM 



relations then give g\ = gz/^/2 = 0.70 ± 0.12 and </ 2 = g^/2 = 1.40 ± 0.24 [30]. However, as discussed in Sec. VI 
the value of "33" extracted directly from r[E* — > A c n] should really be considered as an effective value of the decay 
coupling constant at ttlq = m c , deviating from the static- limit axial coupling by corrections of order Aqcd/ti 
30%. 



VI. CALCULATION OF DECAY WIDTHS 

In this section, we use our lattice QCD results for the axial couplings gi and 93 to calculate various decay widths 
of heavy baryons. At leading order in the chiral expansion, the widths for the strong decays S —¥ T tt are 



T[S->Tx]=$ 



1 



67T/J 



93 



Kj_ 
mn 



^IP I 3 
M S lP7r| ' 



(102) 



where S and T now denote physical si — 1 and s; = heavy-baryon states such as Eb and A&, |p ff | is the magnitude 
of the pion momentum in the S rest frame, 



IPttI 



VWs 



ml][{M s + M T y 



2M.s 



(103) 



and Cf is a flavor factor [HT] , 



Cf = 



1 for 


EW->Ag7r± 


1 for 


E« -► A Q tt , 


1/V2 for 

1 10 f™- 



(104) 



The toq = 00 expression for T can be found for example in |52| . In Eq. ( 102 ), we included the term kj/ijiq to account 
for the first-order corrections for a finite heavy-quark mass. The parameters kj are related to the additional couplings 
in the order-1/mQ HHxPT Lagrangian [S5]. Terms suppressed by (m^/A x ) 2 and (Aqcd/"1q) 2 , which are omitted 
from (102|, lead to small systematic uncertainties in F. 

To determine k x / 2 and k 3 / 2 , we performed fits of experimental data [63] for the widths of the E+ + , E° (J = 1/2) and 
the E* ++ , E*° (J = 3/2) using (102), where we constrained 53 to our lattice QCD result (97) and set itlq = -Mj/^. 



These fits are shown in Fig. H5] and gave the results 



«l/2 
«3/2 



0.55(21) GeV, 
0.47(21) GeV, 



Cov(ki/ 2 , ,g 3 ) = -0.025 GeV, 
Cov(k 3/2 , g 3 ) = -0.025 GeV. 



(105) 



The fit parameters kj are c orrelated with (73, and therefore we also s how the covariances in Eq. (105). The value of 



in Eq. ( 102 ) is plotted as a function of 1/toq in Fig. 16 For tuq 



the sum (73 

are determined dominantly by the experimental input used to fit kj: 

«l/2 



2 M J/i" 



the values of g^ 



rriQ 



33 + 



93 



ifnj/iji 

K 3/2 



■j/4, 



1.059(49), 
= 1.008(46). 



Using the masses of the S* and S c bary ons fr om Ref. 



-> S°7T E+7T 



as shown in Table 



XIII 



we obtain predictions for T[S* 



— > ZjJtt 



Z°TT 4 



(106) 



and 



There, we also show other predictions from the literature, as well as 
upper limits from experiments [71] [75]. Our results for r[S* + ] and r[S*°] are compatible with these limits. 
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FIG. 15. Experimental data for T[Ec — ¥ A c 77*] from Ref. [63], along with fits using Eq. (102 1, for J = 1/2 (solid curve) and 
J = 3/2 (dashed curve). 
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FIG. 16. Value and uncertainty of the quantity (gs + kj/ttiq), which enters in the strong-decay width ( 102 1, as a function of 
the inverse heavy-quark mass ran 1 , for J = 1/2 (solid curve) and J = 3/2 (dashed curve). At m^ 1 = the function is equal 
to 33, which is given by our lattice QCD result (97 1. The vertical lines indicate our choices for the inverse bottom and charm 
quark masses. 



We can also make predictions for the radiative decay S*° — > S° 7, which is forbidden at tree level but can be 
mediated by loops because of flavor-S£/(3) breaking. Using HHxPT, it has been shown that the branching fraction 
of this decay is related to the axial coupling g^ as follows 73 : 



3 „2 



£[Sf ^ ^7] = (1.0 ± 0.3) x 10-^ 



(107) 



Combining this with our lattice QCD result for 52, Eq- (97), and our calculated strong-decay width F[S*° 
S°7T°, S+ 7T-] = 2.78(29) MeV, we obtain 



S[S;°^2° 7 ] = (7±4)xl0- 4 , 
T[S*° -► ~° 7 ] = (2.0 ±1.1) keV. 



(108) 



Next, we discuss the strong decays of bottom baryons. To calculate these widths, we evaluated ( 102 1 for uiq — |Mx 
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In this case the values of <?3 



are determined dominantly by the lattice result (97) for g^: 



H- 



93 



93 



1/2 



= 0.822(87), 



{^ = 0.805(87). 

5 m T 



(109) 



Our calculated widths rp^ > A^tf*] as functions of the X^ — A;, mass difference are shown as the curves in Fig. 17 



Using the experimental values of the baryon masses [131 [33], we obtain the results for T[E 
XIII in agreement with the widths measured by the CDF collaboration [13] . 
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FIG. 17. Widths of the decays T, b * — > Aj, tv as functions of the T. b — Aj, mass difference. The curves (solid: £&, dashed: 
££) and shaded regions show our predictions and their uncertainties. The experimental data points are from CDF [13] . 
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0.51(16) 


2.1 ±1.7 [70] 


1 — 'c 


1.2- 4.1 


1.81 


3.04(37) 


3.18(10) 


2.7(2) 




1.13 


2.44(26) 


< 3.1 (CL=90%) [7TJ 


1 — 'c 


1.2- 4.0 


1.88 


3.12(33) 


3.03(10) 


2.8(2) 




1.08 


2.78(29) 


< 5.5 (CL=90%) [72] 



TABLE XIII. Results in MeV for the total strong-decay widths of charm and bottom baryons. 

In our previous work [5T] we predicted that the widths of the E' b and Sj* are less than 1.1 and 2.8 MeV, respectively. 
Very recently, the CMS collaboration has observed the ££°, finding a width of 2.1 ±1.7 MeV [70fM The mass difference 
to the Sr was measured to be 



M B .» - M = - 



154.41 ± 0.79 MeV. 



(110) 



The E* b ° can decay into S~tt+ (seen by CMS) and into E° b ir°. Taking M s - - M 3 o = 3.1 ± 5.6 ± 1.3 MeV from the 
CDF measurement reported in Ref. [74], we have 

M s ,o - M H o = 157.5 ± 5.8 MeV. (Ill) 



1 Without a spin identification, there is a small possibility that the state observed by CMS is the B.' b ° instead. We do not consider this 
further. 



:-!() 



Using the results ( 110 ) and (111 I, we can update our calculation of the ^° width and find 



b 

T[E* b ° -> S-7T+, -£ tt ] = 0.51 ± 0.16 MeV. (112) 



Given the observed mass difference ( 110 ), and assuming that M^* — M s > « Ms* — M^ b — 21 ± 2 MeV [5S], it is likely 
7/0 _. ct- 



that the decay S 6 ° — >• s b tt + is kinematically forbidden. 



VII. CONCLUSIONS 

The chiral dynamics of mesons and baryons containing a heavy quark is controlled at leading order by three axial 
couplings <7i, 52, and 53. Knowledge of the values of these couplings is an essential ingredient for precision QCD 
calculations in flavor physics. In this paper, we have discussed in detail the first complete lattice QCD determination 
of <7i, 52, and 173. We have extracted the axial couplings by fitting numerical data for matrix elements of the axial 
current using the quark-mass and volume dependence calculated in SU(A\2) heavy-hadron chiral perturbation theory. 
Our final results are 

9l = 0.449 ± 0.047 stat ± 0.019 syst , 

.92 = 0.84 ±0.20 s tat ±0.04 syst , 

g 3 = 0.71 ±0.12 stat ±0.04 syst . (113) 

The systematic uncertainties in ( 113| are very small, because our analysis is based on data at low pion masses, with 



a large volume, and at two different lattice spacings. We have also carefully removed the excited-state contamination 
in the matrix elements by extrapolating the ratios of correlation functions to infinite source-sink separation. 

Previous lattice calculations of heavy-hadron axial couplings had only considered the mesonic coupling g-y. The 
early calculations of g\ did not include dynamical quarks and hence are contaminated by uncontrolled systematic 
errors. The rif = 2 calculations typically used large quark masses and the fits to the quark-mass-dependence were 
performed either linearly in mj or with an incorrect coefficient of the chiral logarithm. Had the correct coefficient 
been used, significantly lower values of g\ would have been obtained in these previous studies. 

For the range of pion masses considered in our work (230 MeV < m w < 350 MeV), the chiral expansion of the 
axial-current matrix elements between heavy-light hadron states is found to be well-behaved. The next-to-leading- 
order contributions are small compared to the leading-order contributions, and NNLO contributions are negligible. 
The rapid convergence of the chiral expansion is also a consequence of the smallness of the static-light axial couplings 



( 113 1. It is interesting to compare the chiral dynamics of hadrons containing a heavy quark with that of light baryons. 
Being particularly light, the interactions of virtual pions (and other pseudo-Goldstone bosons) produce significant 
contributions to many properties of baryons, and generically these effects scale quadratically with the strength with 
which a given baryon sources pions. This, in turn, is determined by the relevant axial coupling, 51,2,3 in the case of 
heavy hadrons, and 17,4 ~ 1.26, |<7ata| ~ 1-6 and 5aa ~ —1.9 in the case of light baryons |75U76j . From the numerical 
values of these couplings, it is apparent that chiral dynamics is more perturbative for heavy-light hadrons than that 
for light baryons. 



Our results for the heavy-light axial couplings, Eq. (113), are significantly smaller than the values one obtains in 



the nonrelativistic quark model, gi — g\ d , g 2 — 2<?^ d and 53 = y/2g^, where g u ^ = 1 is the axial coupling of the 
single-quark transition u — > d. Even if g\ d is set to 0.75, as needed to reproduce the experimental value of the nucleon 



axial charge, the corresponding quark-model values of 51,2,3 are still significantly smaller than the results (113) from 
first-principles lattice QCD. 

We have used our results for g 2 and 53 to calculate strong and radiative decay widths of charm and bottom 
baryons. For the strong decays, we have taken into account the order-1/mQ corrections, which we have constrained 
by combining experimental data for charmed baryon decay rates with our lattice determination of 53. We found that 
the 1/toq corrections are significant (their effect on the amplitudes for Eq — > Aq tt decays is about 40% at mq = m c 
and about 13% at rriQ = nib). As a consequence, the coupling 53 cannot be reliably extracted from experimental 
data for charmed baryon decays alone, and our lattice calculation in the static limit is crucial to calculate the widths 
of bottom baryons. Our results for the widths of the E b baryons are in agreement with recent measurements at 
Fermilab. 

Our determination of the axial couplings can also improve the precision of future lattice QCD calculations of other 
heavy-hadron properties such as masses, decay constants, and form factors, because the axial couplings control the 
dependence of these properties on the light-quark masses. Therefore, the calculation of the axial couplings from first 
principles also has an impact on searches for beyond-the-Standard-Model physics at the LHC and the planned SuperB 
experiment. Importantly, our results include the baryonic couplings 52 and 53. Heavy baryons may offer additional 
opportunities for probing the structure of new physics as a consequence of the different spin quantum numbers. 
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FIG. 18. Summary of all data points for Ri(t), ife(i) and Rs{t). At each value of t/a, results from up to five different values 
of riHYP are shown (from left to right: whyp = 1, 2, 3, 5, 10; points offset horizontally for legibility; in some cases there are no 
results for R3 for the lowest values of hhyp, because the statistical fluctuations were too large to calculate the square root of 
the double ratio). In physical units, the range of the horizontal axis in all plots is from t — 0.336 fm to t — 1.23 fm. 
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Appendix B: Comparison of standard ratio method and summation method 



To extract the effective axiai coupiings from the ratios Ri(t, t 1 ) defined in equations (63), (64), and (65), we defined 
Ri(t) to be the average of Ri(t 1 t') over a symmetric range of t! values around t/2 in a region where there was no 
discernible t' dependence, which essentially amounts to using 



Ri(t, t/2). 



(Bl) 



An alternative approach for extracting g e g is the summation method |19[ I78H81J . In the following, we only consider 
the case of the simple ratios ( 63 ) and ( 64 ) for degenerate spectra. One defines the summed ratio Si (t) by summing 
Ri(t, t') over all values of t', 



Si(t) 



Y,Ri(t,t'). 



(B2) 



For large t, one expects [78j [79 



Si(t) 



(9i 



Icff ' 



(B3) 



with some constant Cj. Thus, the coupling (<?i) c ff can be extracted by taking the derivative |19U81j . 



RT m (t) 



At 



Si(t), 



(B4) 



which is approximated by a finite difference on the lattice. Assuming that there is a nonvanishing off-diagonal matrix 
element of the axial current between the ground-state hadron and an excited state with an energy gap 5 (for our data, 
contamination from off-diagonal matrix elements actually appears to be very small, as discussed in Sees. |IV A| and 
IV B ), one expects that the systematic uncertainties of (JbTJ) and (|B4[) due to this excited state are of order |81J 



Ri(t) 

RT m (t) 



(ffi)eff 

(ft) off 



0(e-i Sit ), 
0(te- 5 ^) 



(B5) 



[see Eq. (71) for the spectral decomposition of Ri(t)]. Thus, the excited-state contamination in i?| ummod (t) decays 
effectively with twice the energy gap relevant for (gi) c s(t), but at the cost of an additional factor of t in front of the 
exponential, which may be important at intermediate values of t. 

Alternatively to taking the derivative as in Eq. ( B4 ), one may fit Si (t) using the linear function ( |B3[ ) with parameters 
Ci and (<7i) c ff- In Fig. 19 we show numerical results for Si(t), along with such fits. In Fig. 20 we compare numerical 
results for the standar d ra tio (Bl), the derivative of the summed ratio (B4), and the results for (gi) s from linear 



fits to Si(t) using Eq. (B3|. For our data, the results from the summation method, especially for the derivative of 
the summed ratio, are seen to suffer from much larger statistical uncertainties than the standard ratio. This was 
also found in Ref. [81] and is not unexpected, because the relative statistical uncertainty in the difference of two 
similarly-sized observables (the discrete derivative used here) is much larger than the relative statistical uncertainty 
in the individual observables. Of course there are correlations which can improve the situation, and we did take these 



into account when calculating ( B4 ) , but because of the way that our lattice calculation was set up (data at successive 
values of t did not always have neighboring source locations), the correlations were not optimal. 

It appears that the systematic errors of the results from the summation method at short t are similar in magnitude 
to the systematic errors of the results from the standard ratio method at the same t, but the deviations from (gi)eS 
have the opposite sign. This shows that valuable information about systematic errors can be obtained by comparing 
both methods. For the present data, our process of extrapolating the results from the standard ratio to infinite t is 
superior because of the much smaller statistical uncertainty. 

Similarly to the work done in Ref. [81] . we also studied models for the three-point and two-point functions with 
excited states. We found that at intermediate values of t, the systematic uncertainties of (Bl ) and (B4 ) were strongly 
dependent on the assumptions made in the model. For some models, the standard ratio showed an adantage while for 
others the summation method showed an advantage, so that again we were not able to draw definitive conclusions. 

Further methods for the calculation of hadron-to-hadron matrix elements are based on the Generalized Eigenvalue 
Problem [5T] and the "Generalized Pencil-of- Function" [52] • These techniques use matrices of correlation functions 
with multiple interpolating fields to reduce the excited-state contamination at finite t. Because we only have data 
from one interpolating field for each hadron, we cannot test these methods here. 
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FIG. 19. Fits to the summed ratios Si(t) and S^i), in the range t/a = 8, 9, 10. The data are for a = 0.112 fm, and a heavy 
quark mass of m„ ^ = 0.04 (close to the physical strange quark mass; the large mass was chosen here for the smaller statistical 
uncertainties) and 71hyp = 3. 
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